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1.  INTRODUCTION 


1.1  Advanced  Gas  Turbine  Requirements 

The  performance  of  the  next  generation  of  military  gas  turbine  engines  will  be  largely 
determined  by  the  temperature  capabilities  of  structural  materials.  This  is  especially  true  in 
the  case  of  the  high  pressure  turbine  where  gas  path  temperatures  already  approach  1800°C 
and  metal  temperatures  approach  1150°C.  Programs  such  as  the  DoD  THPITET  initiative  are 
aimed  at  doubling  the  thrust-to-weight  performance  of  military  engines.  To  accomplish  this 
objective,  design  calculations  indicate  that  gas  path  temperatures  of  about  2100°C  are 
required.  The  same  goal  can  be  accomplished  utilizing  advanced  cooling  schemes  (possible 
with  conductive  materials  such  as  metals  and  intermetallics)  that  maintain  a  gas  path 
temperature  of  2100°C  and  a  metal  temperature  of  1540°C.  These  temperature  requirements 
are  clearly  beyond  the  capabilities  of  materials  currently  available.  In  addition  to  the 
temperature  requirements,  other  critical  materials  properties  cannot  be  compromised.  These 
properties  include  environmental  resistance  (oxidation  and  hot  corrosion),  adequate  low 
temperature  fracture  toughness  or  ductility,  elevated  temperature  strength  (tensile  and  creep), 
and  thermal-mechanical  fatigue  resistance. 

1.2  Materials  Considerations 

The  materials  used  exclusively  in  the  high  pressure  turbine  of  all  state-of-the-art  military 
engines  are  nickel-based  superalloys.  These  materials  are  unique  in  that  they  exhibit 
sufficient  mechanical  properties  and  environmental  resistance  to  be  used  at  temperatures  up  to 
about  85%  of  the  absolute  melting  point.  This  enhanced  capability,  relative  to  other  alloy 
systems,  can  be  attributed  to  the  use  of  coherent  precipitation  strengthening,  alloy 
modifications  with  rhenium  additions  for  enhanced  creep  resistance  and  technologies  such  as 
directional  solidification  for  enhanced  thermal  fatigue  resistance.  Since  the  melting  points  of 
these  alloys  are  about  1350°C,  close  to  that  of  pure  nickel,  significant  further  advancements  in 
temperature  capability  will  not  be  attainable. 

Several  classes  of  materials  can  potentially  meet  the  aggressive  goals  stated  above.  They 
include  ceramics  and  ceramic  matrix  composites  (CMCs),  refractory  metals  with  protective 
coatings,  and  intermetallics  and  intermetallic  matrix  composites  (IMCs).  However,  all  of  these 
materials  will  require  major  improvements  in  one  or  more  critical  properties  to  attain  the 
required  blend  of  environmental  resistance,  low  temperature  ductility  (or  toughness)  and  high 
temperature  strength. 

In  general,  the  ceramics  and  CMCs  suffer  from  a  lack  of  ductility  and  fracture  toughness 
throughout  the  temperature  range.  Additionally,  their  low  thermal  conductivity  would  have  a 
negative  impact  on  the  ability  to  utilize  advanced  cooling  schemes.  They  have  the  advantage 
of  low  densities  and  good  oxidation  resistance.  The  shortcoming  of  available  refractory 
metals  is  poor  oxidation  resistance.  Although  they  can  be  improved  with  protective  coatings 
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(such  as  plasma  sprayed  MoSi2  coatings  or  diffusion  slurry  silicide  coatings  for  niobium 
alloys),  failure  of  these  coatings  results  in  rapid  deterioration  of  the  substrate  due  to  either 
internal  oxidation  (resulting  in  alloy  embrittlement)  or  non-protective  oxidation.  Refractory 
metals  have  the  advantages  of  having  a  mature  metallurgical  knowledge  and  experience  base, 
good  thermal  conductivities,  adequate  ductility  throughout  the  temperature  range  of  interest, 
and  can  be  strengthened  by  traditional  solid  solution  and  dispersion  hardening  techniques'. 
Intermetallics  generally  exhibit  insufficient  ductility  at  low  temperatures  (below  the  alloy 
ductile-to-brittle  transition  temperature  [DBTT])  and  inadequate  strength  at  elevated 
temperatures.  They  have  good  thermal  conductivities  such  that  advanced  cooling  schemes  can 
be  used  and  they  have  lower  densities  than  the  refractory  metals  (but  generally  higher  than  the 
ceramic  materials).  Several  of  the  intermetallic  systems  such  as  NiAl  and  MoSi2  exhibit 
excellent  oxidation  resistances.  Of  all  the  intermetallics,  MoSi2  has  the  best  combination  of  a 
high  melting  point,  high  temperature  strength,  and  oxidation  resistance. 

1.3  MoSii  Background 

The  properties  of  MoSi2  have  been  investigated  for  potential  applications  as  a  high 
temperature  structural  material,  furnace  heating  elements,  protective  coatings  for  refractory 
metals,  and  thin  film  diffusion  barriers  at  interconnects  in  semiconductor  devices.  The 
properties  of  interest  for  potential  application  in  advanced  gas  turbine  engines  are  (1)  elevated 
temperature  oxidation  resistance  (2)  intermediate  temperature  oxidation  resistance  (3)  high 
temperature  strength  (including  creep)  (4)  brittle-to-ductile  transition  and  (5)  low  temperature 
fracture  toughness. 

The  acceptance  of  MoSi2  as  a  high  temperature  material  is  primarily  attributable  to  its 
excellent  oxidation  resistance.  It  is  the  material  of  choice  as  heating  elements  for  air  furnace 
applications  to  temperatures  slightly  in  excess  of  1800°C.  This  remarkable  protection  is  due 
to  the  ability  of  the  material  to  form  a  pure  Si02  scale.  Excellent  oxidation  resistance  is 
attained  because  of  three  factors.  First,  the  free  energy  of  formation  for  Si02  is  significantly 
lower  than  for  M0O3  and  favors  Si02  formation.  Second,  M0O3  has  a  very  high  vapor 
pressure  and  at  temperatures  above  about  800°C,  will  evaporate  and  leave  behind  a  pure  Si02 
scale.  Third,  the  growth  rate  of  the  Si02  scale  is  very  slow.  The  parabolic  rate  constants  (a 
measure  of  scale  growth  rate)  shown  in  Fig.  1  indicate  a  Si02  scale  grows  more  slowly  than  a 
AI2O3  scale  at  temperatures  above  about  1500°C.  The  oxidation  resistance  of  MoSi2  at 
intermediate  temperatures  is  of  some  concern.  At  temperatures  ranging  from  500  to  600°C  ,  a 
phenomenon  called  "pesting"  has  been  observed.  This  phrase  was  first  referenced  by  Fitzer^ 
and  is  generally  defined  as  the  disintegration  of  a  material  after  cyclic  exposure  to 
intermediate  temperatures.  More  recent  work  on  MoSi2  by  most  investigators  have  not 
indicated  a  pesting  problem.  This  suggests  that  the  problem  is  associated  with  the  processing 
and/or  purity  of  early  material.  This  was  confirmed  by  Meier^,  who  observed  pesting  in  arc- 
melted  and  porous  P/M  MoSi2  after  exposures  at  500°C. 
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The  mechanical  properties  of  monolithic  MoSia  indicate  the  material,  in  pure  form,  exhibits 
inadequate  ductility  or  fracture  toughness  at  low  temperatures  and  insufficient  strength  at 
elevated  temperatures.  The  properties  can  be  characterized  by  three  different  performance 
regions'^.  At  low  temperatures  (below  the  DBTT  of  about  1000°C),  MoSi2  is  strong  and 
brittle.  At  intermediate  temperatures  (above  the  1000°C  to  about  1250°C)  MoSi2  is  strong 
and  ductile.  At  elevated  temperatures  (above  about  1250°C),  MoSi2  is  ductile  but  exhibits 
low  strength.  The  ultimate  strength  of  pure  MoSi2  has  been  determined  by  various 
investigators  and  is  constant  at  280-350  MPa  from  low  temperatures  to  about  1250°C  as 
shown  in  Fig.  2.  Above  this  temperature,  the  strength  decreases  rapidly  (due  to  thermally 
activated  processes).  This  strength  level  is  lower  than  that  of  cobalt  and  nickel-based 
superalloys  at  low  temperatures  where  strengths  of  about  500  to  1000  MPa  is  required  for 
highly  stressed  components  such  as  turbine  blades  and  vanes.  At  temperatures  above  1 100°C 
when  the  strength  of  the  cobalt  and  nickel  alloys  decreases  rapidly  as  we  approach  their 
melting  points,  the  MoSi2  exhibits  higher  strengths.  Similarly,  MoSi2  exhibits  greater  creep 
resistance  above  1 100°C  as  shown  by  the  results  in  Fig.  3.  The  fracture  toughness  of  MoSi2  is 
only  about  2  MPa  Vm  which  is  typical  of  most  intermetallic  and  ceramic  materials.  These 
mechanical  property  results  certainly  indicate  a  need  to  improve  the  low  temperature  fracture 
toughness  and  elevated  temperature  mechanical  properties  of  MoSi2. 

Only  limited  investigations  have  been  found  which  dealt  with  improving  the  mechanical 
properties  of  MoSi2  through  alloying.  Researchers^  have  found  that  the  addition  of 
germanium  to  MoSi2  improves  the  cyclic  oxidation  resistance  of  MoSi2  and  is  especially 
effective  in  minimizing  pesting.  Altering  the  viscosity  and  coefficient  of  thermal  expansion  of 
the  Si02  scale  results  from  germanium  additions.  This  minimizes  thermal  stresses  between 
the  alloy  and  the  oxide  scale,  reducing  scale  cracking  and  spallation.  To  improve  the 
performance  of  heating  elements,  Kanthal  Corp.  alloys  MoSi2  with  significant  levels  of 
tungsten.  An  examination  of  the  Mo-W-Si  ternary  phase  diagram  indicates  that  alloying  with 
about  50  atom  percent  tungsten  for  molybdenum  would  improve  the  melting  point  and  the 
high  temperature  capabilities  of  MoSi2  by  about  70°C  (thus  extending  the  use  temperature  to 
about  1900°C).  In  work  performed  at  UTRC  that  was  internally  funded  and  subsequently 
sponsored  by  WRDC^,  rhenium  modifications  were  found  to  significantly  reduce  the 
minimum  creep  rate  of  MoSi2.  The  results  in  Fig.  4  indicate  the  minimum  creep  rate  for 
rhenium  modified  MoSi2  exhibits  one  to  two  orders  of  magnitude  lower  creep  rate  than 
unalloyed  MoSi2  at  1000°C  and  at  1400°C.  In  this  same  study,  tungsten  was  found  to  have 
little  effect  on  elevated  temperature  creep  strength.  An  understanding  of  these  compositional 
effects  has  not  been  established. 

Intermetallic  matrix  composites  based  on  MoSi2  are  being  investigated  to  meet  the 
requirements  of  the  next  generation  gas  turbine  engines.  The  objectives  of  these  investigators 
are  the  same  as  ours,  that  is,  to  improve  the  low  temperature  fracture  toughness  and  elevated 
temperature  strength  of  MoSi2.  Most  of  the  work  is  currently  centered  on  the  fabrication  of 
MoSi2  based  intermetallics  and  on  the  effect  and  compatibility  of  different  reinforcing  fibers. 
The  SiC  whisker/MoSi2  matrix  composite  system  has  attracted  the  most  interest.  This  is 
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attributed  to  the  availability  of  the  reinforcements  and  to  the  oxidation  resistance  and  thermal 
stability  of  the  SiC  with  respect  to  MoSi2.  The  results  of  the  various  studies  indicate  only 
modest  improvements  in  strength  and  fracture  toughness  were  attained.  For  example,  the 
fracture  toughness  of  MoSi2  was  improved^  from  2  to  only  about  5  MPa  Vm.  Similarly,  only  a 
marginal  improvement  in  elevated  temperature  strength  was  observed.  The  greatest  increases 
in  fracture  toughness  have  been  obtained  utilizing  refractory  metal  fibers.  For  example, 
fracture  toughnesses  of  up  to  39  MPa  Vm  have  been  achieved*.  Unfortunately,  refractory 
metal  fibers  are  not  thermodynamically  stable  with  respect  to  MoSi2  and  either  dissolve  (and 
alloy  with  MoSi2  to  form  a  lower  silicide)  or  react  to  form  a  more  stable  silicide  such  as 
NbSi2.  Additionally,  any  exposed  refractory  metal  fibers  will  compromise  the  oxidation 
resistance  of  MoSi2.  The  current  work  on  alloy  matrix  composites  have  been  concentrated  in 
the  (W-Mo)Si2  matrix/SiC  fiber  system.  As  expected,  only  modest  improvements  in  strength 
and  fracture  toughness  have  been  achieved  using  this  alloy  matrix. 

These  results  indicate  MoS^  offers  excellent  oxidation  resistance  along  with  the  potential  for 
high  temperature  strengthening  through  compositing  and  matrix  alloying.  With  respect  to 
alloying,  little  has  been  done  to  investigate  compositional  effects  on  the  mechanical  properties 
of  MoSi2.  This  program  is  aimed  at  improving  and  understanding  the  mechanical  properties 
of  MoSi2  through  compositional  modifications  and  towards  the  understanding  and 
development  of  materials  having  the  balanced  properties  needed  in  the  next  generation  gas 
turbine  engines. 


2.  ALLOY  SELECTION 


2.1  Approach 

In  this  program,  the  compositional  effects  on  the  mechanical  properties  of  MoSia  are  being 
investigated.  The  objective  of  this  activity  is  to  predict,  study  and  reconcile  alloying  effects 
on  the  elevated  temperature  strength  and  low  temperature  fracture  toughness.  The  approach 
used  is  to  select  appropriate  alloying  elements  based  on  state-of-the-art  alloy  design  theory. 
After  suitable  alloys  have  been  identified,  materials  fabrication  for  evaluation  is  accomplished 
by  appropriate  arc-melting,  powder  metallurgy  (including  vacuum  hot  pressing  or  hot  isostatic 
pressing)  and  plasma  spraying  techniques  which  have  already  been  developed  for  preparation 
of  intermetallic  alloys.  The  phase  stability,  structure,  and  mechanical  properties  are  being 
measured,  correlated  with  predicted  effects  and  compared  to  the  baseline  MoSi2  properties. 

The  significance  of  this  work  is  not  restricted  to  achieving  a  fundamental  understanding  of 
alloying  effects  on  MoSia,  but  will  also  be  realized  when  used  in  conjunction  with  the  best 
compositing  techniques  being  evaluated  by  other  investigators.  It  is  viewed  that  this  program 
would  complement  the  composite  effort  on  improving  low  temperature  fracture  toughness  and 
elevated  temperature  strength.  Both  of  these  efforts  are  necessary  to  achieve  the  goal  of 
developing  suitable  materials  for  the  next  generation  gas  turbine  engines. 

2.2  Alloy  Theory 

The  selection  of  alloying  elements  for  evaluation  is  based  on  the  best  available  phase 
equilibria,  alloy  theory,  simulation  or  modeling  capabilities,  and  property  requirements.  Phase 
equilibria  and  crystal  structures  in  the  Mo-Si-X  ternary  system  is  of  primary  importance. 
Elemental  additions  that  compromise  the  melting  point  (by  forming  lower  melting  silicides) 
were  not  selected.  Since  most  structural  materials  are  usable  to  about  O.S-O.bTj^,  any 
significant  reduction  in  the  melting  point  of  MoSia  (2020°C)  is  undesirable.  Based  on  this 
criterion,  the  best  elemental  candidates  are  tungsten,  rhenium,  tantalum,  niobium,  germanium 
and  carbon.  The  solubility  limits  of  ternary  alloying  additions  into  MoSia  are  not  well 
documented  and  while  a  number  of  these  diagrams  are  published,  there  are  numerous 
inconsistencies  and  errors.  For  example,  the  ternary  Cr-Mo-Si  phase  diagram’  shown  in  Fig. 
5  indicates  a  large  MoSij  phase  field  with  the  Cllb  crystal  structure  but  internal  work  at 
UTRC  indicates  this  solubility  range  only  extends  to  about  2  weight  percent  chromium.  More 
detailed  ternary  phase  diagrams  are  being  made  of  each  alloy  system  studied. 

The  crystal  stmcture  of  a  material  has  a  significant  effect  on  the  mechanical  properties. 
Differences  in  the  crystal  structure  will  have  significant  impact  on  the  number  of  slip  systems 
available  and  the  operative  slip  vectors.  In  addition  to  solid  solution  alloying,  additions  will 
be  selected  to  stabilize  specific  crystal  structures  in  the  MoSi2  system.  The  equilibrium 
structure  for  MoSi2  is  the  CII5  tetragonal  phase.  There  is  currently  a  controversy  as  to  the 
existence  of  a  stable  C40  hexagonal  stmcture  in  the  MoSi2  system.  This  stmcture  is  either  a 
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metastable  phase*”  or  a  high  temperature  phase**  (occurring  above  1900°C).  Interstitial 
elements  such  as  carbon  are  believed  to  stabilize  the  C40  phase  and  research  is  being 
conducted  in  this  area*^. 

Engel-Brewer  diagrams*^  provide  a  perspective  on  the  nature  of  crystal  structure  similarities 
which  are  not  evident  from  crystallographic  data;  such  data  consists  primarily  of  symmetry 
elements  and  overlooks  relative  differences  in  structure.  The  strength  of  these  diagrams  over 
other  mapping  schemes  is  that  the  solubility  limits  of  compounds  and  solid  solutions  are 
readily  apparent  as  a  function  of  ternary  alloying.  These  maps  are  simply  composites  of 
binary  phase  diagrams  taken  at  one  temperature.  An  examination  of  such  a  map  with  silicon 
indicates  the  Cllt,  structure  is  stable  for  Mo,  W,  and  Re.  One  notes  that  disilicides  with 
lower  electron/atom  ratios  than  MoSii  (Cllj,)  have  the  C40  structure.  The  C40  structure  is 
simply  a  shear  vector  removed  from  the  Clli,  structure*^.  Other  closely  related  structures 
such  as  C54  and  C49  have  also  been  identified.  While  both  Engel-Brewer  and  Pettifor*'*  maps 
would  indicate  that  the  C40,  C49,  and  C54  crystal  structures  are  closely  related  to  that  of 
Cll|j,  the  actual  geometric  relationship  has  only  been  obtained  for  the  C40  structure.  These 
relationships  have  not  been  explored  in  great  detail  and  simply  detailing  the  atomic  motions 
required  for  the  transformations  will  be  a  significant  task. 

Having  isolated  the  structural  similarities,  partial  destabilization  of  the  Cllb  structure  can 
lead  to  coherent  precipitation  or  enhanced  ternary  solubility  which  could  lead  to  greater  creep 
resistance  while  a  lowered  ordering  parameter  may  result  in  a  lowered  DBTT  and  enhanced 
low  temperature  ductility.  An  example  of  the  destabilization  of  the  Cl  Ij,  structure  leading  to 
enhanced  ternary  solubility  occurs  in  the  case  of  chromium  additions  to  MoSi2.  CrSi2  has  the 
C40  structure  and  is  expected  to  partially  destabilize  CII5  towards  the  C40  structure. 
Examination  of  this  phase  diagram  would  indicate  that  the  MoSi2  phase  field  can  be  expanded 
considerably  in  ternary  space.  This  is  attributed  to  chromium  substitution  for  either 
molybdenum  or  silicon  in  MoSi2  while  molybdenum  substitutes  only  for  chromium  in  CrSi2. 
The  effect  on  strength  of  MoSi2  through  chromium  additions  has  been  recently  documented*^. 
It  is  expected  that  other  such  alloying  approaches  can  be  used  to  strengthen  MoSi2. 

Similarly,  titanium  is  expected  to  stabilize  the  C54  structure  (TiSi2)  while  zirconium  or 
hafnium  is  expected  to  stabilize  the  C49  structure  (ZrSi2  and  HfSi2),  and  their  additions  are 
being  considered. 

With  respect  to  major  alloying  modifications,  the  two  methods  most  cited  in  the  literature*^’*’ 
as  being  responsible  for  altering  the  crystal  structure  and  properties  of  intermetallic 
compounds  are  (1)  electron-to-atom  ratios  (e/a)  of  the  compounds*®’**  and  (2)  relative  atomic 
radii  ratios  (Ra/Rb)  of  the  compounds*’’*^.  These  are  accomplished  through  substitutional 
alloying.  The  e/a  ratio  modifications  are  accomplished  by  substituting  elements  into  MoSi2 
which  are  expected  to  maintain  single  phase  compounds  but  which  contain  either  an  excess  or 
deficit  of  electrons  with  respect  to  the  substituted  element.  Determination  of  electron 
densities  per  atom  is  done  by  counting  the  average  electrons  per  atom  in  excess  of  the  inert 
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gas  core  electrons.  This  has  been  shown  to  be  an  adequate  measure  of  electron  density^®.  For 
example,  substitution  of  rhenium  for  25  atomic  percent  molybdenum  would  increase  the  e/a 
from  4.67  to  4.75.  Similarly,  substitution  of  niobium  for  molybdenum  would  decrease  the  e/a 
from  4.67  to  4.58. 

Atomic  radii  ratio  modifications  are  accomplished  in  a  similar  way  but  with  the  e/a  ratio 
constant.  Since  the  atomic  radii  in  intermetallic  compounds  are  not  fully  understood,  the  hard 
sphere  approximation  is  being  used.  More  suitable  radii  will  also  be  sought  such  as  ionic  radii 
or  some  fraction  of  atomic  radii  which  may  yield  results  that  can  be  correlated  with  the 
experimental  findings.  As  an  example,  an  alloy  modification  to  MoSia  to  change  Ra/Rb 
would  be  to  substitute  germanium  which  has  an  atomic  radius  of  0.139nm  for  silicon  which  as 
an  atomic  radius  of  0. 1 17nm. 

Based  on  the  approaches  described,  six  elemental  modifications  were  selected  for  study  in  the 
initial  phase  of  this  program.  The  modifications  were  selected  based  on  available  phase 
diagrams,  melting  temperatures,  crystal  structures,  e/a  ratios  and  atomic  radii. 

The  elements  selected  are  B,  Al,  Ge,  Re,  Nb,  and  Hf.  Boron  was  selected  in  an  attempt  to 
stabilize  the  C32  structure  as  well  as  evaluating  the  effect  of  smaller  atomic  radii  (boron  has 
an  atomic  radius  of  0.097nm  while  silicon  has  an  atomic  radius  of  0.1 17nm).  Aluminum  was 
selected  because  of  the  large  solubility  range  of  Mo(Si,Al)2  which  has  the  C40  structure. 
Germanium  was  selected  because  it  has  complete  solubility  in  MoSii  with  the  Cl  Ij,  structure 
when  substituted  for  silicon.  Gemanium  also  has  a  significantly  larger  atomic  radius 
(0.139nm)  than  silicon  (0.117nm).  Re  was  selected  because  of  its  complete  solubility  in 
MoSi2  with  the  Cll|j  structure.  Additionally,  the  melting  temperature  and  e/a  ratio  are 
expected  to  increase  with  rhenium  substitutions  for  molybdenum.  Niobium  was  selected 
because  a  high  solubility  of  Nb  in  MoSi2  (Cl  lb)  and  Mo  in  NbSi2  (C40)  are  expected, 
additionally,  the  e/a  ratio  is  reduced  with  niobium  substitutions  for  molybdenum.  Hafnium 
modifications  were  selected  to  attempt  to  stabilitze  the  C49  crystal  structure  as  well  as  further 
reducing  the  e/a  ratio  with  hafnium  substitutions  for  molybdenum. 

2.3  Atomic  Modeling 

Materials  properties  controlled  at  the  atomic  level  can  in  principle  be  calculated  from  quantum 
mechanical  methods.  These  phenomena  include  the  ground  and  excited  crystal  structures 
along  with  their  lattice  parameters  and  enthalpies  of  transformation,  coefficients  of  thermal 
expansion,  single  crystal  elastic  constants,  PiereTs  barrier,  dislocation  core  energies,  and 
barriers  to  diffusion.  While  there  are  no  limitations,  in  principle,  to  the  accuracy  with  which 
these  calculations  can  be  performed,  there  are  a  number  of  practical  considerations  which 
have  complicated  the  uses  of  these  techniques  as  metallurgical  tools. 
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Among  the  most  significant  practical  limitations  are  those  concerned  with  the  numerical 
solution  of  the  Schroedinger  Wave  Equation  for  a  large  number  of  symmetry-unique  atoms. 
This  limitation  has  motivated  a  variety  of  approximate  methods  which  are  now  employed  to 
generate  solutions  to  this  equation.  As  the  speed  of  digital  computers  has  continued  to 
increase,  the  severity  of  the  approximations  required  has  diminished.  Today  the  ready 
availability  of  quantum  mechanical  software,  coupled  with  high  speed  workstations  makes  it 
possible  that  metallurgical  research  can  prosper  from  the  use  of  quantum  mechanical  methods. 
There  is  also  every  reason  to  believe  that  as  the  number  of  investigators  using  these  methods 
increases,  new  capabilities  will  be  provided  at  a  rapid  pace.  Some  of  the  underlying 
principles  of  quantum  mechanical  methods  will  be  summarized  here  in  an  attempt  to  give  the 
average  metallurgist  a  firmer  foundation  from  which  to  appraise  the  usefulness  of  these 
methods  in  his  or  her  research  program. 

To  begin  with,  all  quantum  mechanical  methods  are  involved  in  finding  the  solutions  to  either 
the  Schroedinger  or  Dirac  equation.  Of  the  two,  the  Dirac  equation  is  far  more  difficult  to 
solve,  but  its  solutions  include  relativistic  effects  experienced  by  electrons  as  they  move  about 
the  nucleus.  For  the  most  part,  these  effects  are  very  small,  but  become  appreciable  as  the 
positive  charge  of  the  nucleus  increases.  These  effects  are  considered  to  be  significant  for 
atoms  heavier  than  Hf.  While  the  solutions  to  the  Schroedinger  equation  do  not  include  these 
relativistic  effects  directly,  corrections  can  be  made  to  include  many  of  these  effects  indirectly. 
The  methods,  which  solve  the  Schroedinger  equation  with  varying  degrees  of  relativistic 
corrections,  are  the  most  common;  it  is  those  solutions  and  methods  which  will  be  discussed 
in  this  report. 

The  solution  to  any  quantum  mechanical  equation,  including  the  Schroedinger  equation,  gives 
two  results.  The  first  of  these  is  an  electronic  energy,  which  is  equivalent  to  the 
thermodynamic  energy,  i.e.  the  enthalpy,  of  the  system  at  zero  degrees  Kelvin.  This  energy  is 
generally  orders  of  magnitude  larger  than  vibrational  and  rotational  energies,  which  become 
appreciable  at  high  temperatures,  and  is  itself  nearly  constant  as  a  function  of  temperature. 
The  other  result  of  the  solution  to  the  Schroedinger  equation  is  a  wavefunction,  which 
provides  information  regarding  the  electron  distribution  in  the  system. 

As  the  energy  resulting  from  a  QM  calculation  represents  an  enthalpy,  the  changes  in  enthalpy 
associated  with  various  structural  and  chemical  changes  can  in  principle  be  determined 
through  quantum  mechanical  methods.  Combined  with  statistical  approaches,  one  can  then 
deduce  a  change  in  free  energy  and  consequently  all  of  the  thermodynamic  properties  of  a 
system  of  interest.  While  the  utility  of  the  quantum  mechanically  calculated  energy  should  be 
obvious,  the  practical  uses  of  the  charge  density  or  wavefunction  is  more  abstract.  Often  this 
information  has  been  utilized  to  provide  insight  about  the  nature  of  the  bonding  present  in  a 
material  or  molecule.  It  should  be  emphasized,  however,  that  the  chemical  bond  is  an 
abstraction  for  which  there  is  no  physical  analog. 


8 


The  single  rigorous  relationship  between  the  charge  density  and  materials  properties  makes 
use  of  the  electrostatic  theorem.  This  theorem  states  that  the  force  exerted  on  a  nucleus  from 
a  given  charge  distribution  can  be  determined  from  classical  electrostatics.  So  once  a  charge 
density  of  purely  quantum  mechanical  origin  has  been  determined,  it  is  possible  to  analyze  the 
force  on  the  nuclei  of  that  system  from  classical  electrostatics.  It  is  from  this  theory  that  much 
of  the  concept  of  bonding  has  been  developed.  In  a  molecular  system  where  there  is  a  buildup 
of  charge  density  between  atoms,  the  atoms  are  said  to  be  bound,  for,  assuming  that  the  charge 
density  is  rigid,  the  local  density  will  act  to  hold  the  nuclei  together  as  one  tries  to  separate 
them. 

Writing  down  the  Schroedinger  equation  for  a  specific  molecule  or  solid  involves  determining 
the  potential  acting  on  each  electron.  There  are  two  parts  to  this  potential:  i)  the  attraction 
between  the  electron  and  each  nuclei,  and  ii)  the  repulsion  between  the  electron  and  all  others. 
This  latter  potential,  electron-electron  repulsion,  gives  rise  to  the  phenomenon  of  electron 
correlation,  whereby  electrons  will  move  in  a  fashion  so  as  to  stay  out  of  each  other’s  way. 
The  electron-electron  repulsion  also  guarantees  that  for  all  systems  with  more  than  two 
electrons,  there  is  no  solution  to  the  Schroedinger  equation  which  can  be  represented  as  a 
simple  analytical  function.  While  it  is,  in  principle,  possible  to  calculate  a  wave  function  and 
energy  for  any  system  exactly,  the  solutions  will  always  take  the  form  of  an  infinite  series.  In 
practice,  therefore,  one  must  truncate  these  series  after  some  number  of  terms,  giving  rise  to 
an  error  in  both  the  energy  and  wavefunction.  There  are  two  approaches  to  minimizing  this 
error.  One,  is  to  find  rapidly  converging  basis  sets  in  which  to  expand  the  solutions,  while  the 
other  is  to  use  basis  sets  which  provide  some  mathematical  convenience  allowing  one  to 
handle  more  terms  in  the  expansion. 

One  of  the  more  rapidly  converging  basis  sets  is  that  made  up  of  the  so  called  spherical 
harmonics.  These  functions  are  the  familiar  s,  p,  d,  and  f  orbitals,  the  shapes  of  which  are 
familiar  to  all  freshmen  chemistry  students.  These  are,  in  fact,  the  angular  part  of  the 
solutions  to  the  central  force  problem,  i.e.  that  part  of  the  atomic  Schroedinger  equation  which 
does  not  involve  electron-electron  repulsion.  This  being  the  case,  it  is  no  wonder  that  these 
functions  are  rapidly  convergent,  at  least  when  electron  correlation  is  ignored. 

When  electron  correlation  is  included,  many  more  spherical  harmonics  must  be  included  in 
the  expansion  of  the  wavefunction.  Including  these  terms  is  referred  to  as  including 
configuration  interactions  or  Cl  terms.  The  larger  the  basis  set  in  which  the  solution  is 
expanded  the  more  computationally  intense  the  problem.  Including  Cl  interactions  rapidly 
increases  the  complexity  of  the  problem.  As  such  this  approach  is  used  primarily  for  small 
molecules  where  a  high  degree  of  accuracy  is  required.  These  approaches  are  referred  to 
collectively  as  Hartree-Fock  (HF)  methods,  and  are  distinguished  by  the  fact  that  one  can 
systematically  improve  the  accuracy  of  the  calculation  by  increasing  the  size  of  the  basis  set. 
In  principle,  HF  methods  can  produce  an  exact  answer,  though  in  application  there  is  always 
some  error  introduced  by  the  limited  basis  set. 
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Another  approach  to  solving  the  Schroedinger  equation  which  can,  in  theory,  produce  an  exact 
answer  is  referred  to  as  a  local-density-functional  (LDF)  approach.  This  approach  make  use 
of  a  theorem  proved  by  Kohn  and  Sham  showing  that  the  energy  of  a  molecular  or  solid  state 
system  is  a  functional  of  the  electron  density.  The  Schroedinger  equation  cast  in  the  form  of  a 
density  functional  equation  is  referred  to  as  the  Kohn-Sham  equation.  The  disadvantage  in 
using  LDF  methods,  is  that  the  exact  functional  form  by  which  the  energy  depends  on  the 
electron  density  is  unknown.  As  a  consequence,  a  number  of  functional  forms  have  been 
suggested  and  are  constantly  being  refined  to  produce  answers  which  are  consistent  with 
experiment  and  other  theories.  The  advantage  of  LDF  methods  is  that  the  results  are  not  as 
sensitive  to  basis  set  size,  with  electron  correlation  and  exchange  included  in  the  functional. 
As  a  result,  smaller  basis  sets  are  required  and  larger  systems  can  be  realistically  treated.  Just 
as  with  HF  methods,  LDF  methods  can  theoretically  provide  exact  solutions  while  in  practice, 
some  error  is  introduced:  not  as  much  by  a  limited  basis  set,  but  by  the  form  of  the  density 
functional.  For  the  most  part  LDF  techniques  have  been  those  most  commonly  used  to  study 
solid  state  phenomena,  while  HF  approaches  have  been  preferred  in  the  study  of  molecular 
properties.  As  LDF  methods  are  not  hampered  by  needing  to  use  extremely  large  basis  sets  to 
treat  exchange  and  correlation,  novel  methods  could  be  used  in  solving  Kohn-Sham  equation 
without  introducing  excessive  additional  errors.  This  was  particularly  true  in  the  60  s  and  70  s 
where  even  the  super-computers  of  the  day  were  limited  in  terms  of  the  size  of  the  problem 
which  could  be  handled,  even  when  employing  the  LDF  approximations. 

One  of  the  novel  methods  most  deserving  of  comment  is  that  involving  the  partitioning  of  the 
molecular  or  solid  state  potential  into  a  spherical  part,  called  the  muffin  tin,  and  a  constant 
part  located  between  atoms  in  the  interstitial  region.  For  the  time,  this  was  an  excellent 
approximation,  since  to  an  electron  near  an  atomic  nucleus  “sees”  primarily  the  central  force 
from  that  nucleus,  while  in  the  interstitial  region  the  potential  is  slowly  varying  and  can  be 
treated  as  constant.  For  a  constant  spherical  potential,  the  solutions  to  the  Schroedinger 
equation  are  in  the  form  of  a  radial  part  (a  linear  combination  of  spherical  harmonics)  near  the 
nucleus  and,  in  the  interstitial  region,  the  solutions  are  in  the  form  of  plane  waves.  By 
requiring  that  the  solutions  between  the  two  regions  be  smooth  and  continuous,  one  generates 
an  acceptable  solution  to  the  molecular  or  solid  state  problem.  These  approaches  are 
collectively  referred  to  as  muffin  tin  techniques  and  assume  some  shape  for  the  potential. 
Today  it  is  no  longer  necessary  to  assume  a  shape  of  the  molecular  or  solid  state  potential  and 
methods  which  do  not  assume  a  potential  shape  are  called  full  potential  techniques.  Even 
though  these  methods  make  no  approximations  to  the  molecular  potential  they  still  often 
partition  space  into  an  atomic  region  and  an  interstitial  region,  separately,  solving  the 
Schroedinger  equation  in  each  region  and  then  matching  solutions  across  the  region 
boundaries  to  produce  the  correct  full  solution. 

All  current  full  potential  LDF  codes  are  capable  of  producing  results  with  similar  accuracy. 
Using  these  full  potential  techniques  along  with  the  newest  local-density-functionals,  current 
calculations  are  capable  of  determining  equilibrium  lattice  constants  within  about  2%,  elastic 
constants  to  better  than  15%;  all  of  this  can  be  done  for  unit  cells  containing  as  many  as  15 
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atoms.  What  distinguishes  one  full  potential  approach  from  another  is  the  basis  set  in  which 
the  solution  is  expanded.  For  example,  plane  wave  expansions  offer  a  number  of  attractive 
mathematical  properties;  consequently,  there  are  a  number  of  methods  which  make  use  of 
plane  wave  expansions.  One  of  the  best  known  of  these  is  the  FLAPW  method,  standing  for 
full  potential  linear  augmented  plane  wave  method.  Another  method  is  the  FLASTO  method, 
which  is  a  full  potential  linear  augmented  Slater  type  orbital  approach.  In  this  approach  Slater 
type  orbitals  are  used  to  expand  solutions  to  the  Schroedinger  equations. 

Other  letters  of  the  alphabet  soup  used  in  referring  to  electronic  structure  methods  relay 
information  regarding  how  a  problem  is  solved  rather  than  the  basis  set  in  which  it  is  solved. 
For  example,  KKR  methods  are  Green’s  function  techniques  which  solve  the  Schroedinger 
equation  as  a  scattering  problem.  Many  of  these  techniques  use  the  muffin  tin  approximation 
e.g.  the  LKKR  technique. 

Another  advantage  of  full  potential  approaches,  in  addition  to  their  superior  accuracy  when 
determining  energies,  is  that  of  producing  smooth  charge  densities  and  wavefunctions.  Unlike 
muffin  tin  approaches,  the  cusps  and  sharp  comers  where  the  wavefunctions  cross  the  muffin 
tin  boundaries  are  not  seen.  For  those  interested  in  establishing  correlations  between 
variations  in  charge  density  and  properties,  full  potential  techniques  should  be  used.  All  full 
potential  techniques  provide  the  necessary  accuracy  to  determine  a  number  of  quantities  of 
interest  to  metallurgists,  what  has  complicated  their  use  is  that  most  really  interesting 
metallurgical  parameters:  defects,  diffusion  barriers,  ternary  and  quaternary  effects,  etc.  can 
only  be  modeled  as  large  supercells,  or  clusters,  with  very  many  symmetry  unique  atoms  per 
unit  cell  or  in  the  cluster.  These  very  large  systems  offer  their  own  complications,  and 
research  has  proceeded  along  two  courses. 

The  first  of  these  hcis  relied  on  the  increasing  speed  of  computers.  Using  the  same  algorithms 
that  15  years  ago  allowed  for  full  potential  calculations  on  three  symmetry  unique  atoms  per 
unit  cell,  calculations  can  now  be  performed  on  unit  cells  containing  ten  or  fifteen  symmetry 
unique  atoms  and  require  approximately  the  same  amount  of  CPU  time. 

Combing  increased  computer  speed  with  improved  algorithms  now  allows  the  modeling  of 
supercells  with  up  to  50  atoms  per  cell.  The  most  important  improvement  in  computational 
approach  comes  from  a  set  of  methods  known  as  linear  methods.  This  is  the  L  in  FLAPW, 
LASTO,  or  LMTO.  These  methods  are  all  accurate  and  fast.  Still  these  approaches  are  not 
capable  of  modeling,  for  example,  a  dislocation  core,  which  might  require  a  1000  atoms  per 
unit  cell.  However,  given  the  speed  of  computer  processors  is  doubling  once  a  year  and  multi¬ 
processor  machines  are  more  conunon-place,  in  time  such  large  calculations  should  be 
practical.  Therefore,  researchers  following  this  course  have  been  involved  in  increasing  the 
speed  of  their  computations,  without  sacrificing  accuracy,  anticipating  that  in  the  not  to  distant 
future  vary  large  and  meaningful  problems  can  be  handled  with  the  techniques  which  can  in 
principle  provide  exact  solutions. 
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There  are  also  some  specialty  codes  which  were  designed  to  treat  a  particular  class  of 
phenomena.  The  LKKR  is  an  example.  This  approaches  allows  one  to  treat  isolated  coherent 
interfaces  and  has  proven  useful  in  the  determining  boundary  energies  such  as  SF’s,  APB’s, 
and  grain  boundaries.  This  technique  uses  the  muffin  tin  approximation  and  therefore  is  best 
suited  to  modeling  boundaries  with  a  minimum  of  free  volume. 

The  second  course  of  research  has  involved  adopting  approximate  methods  which  allows  one 
to  model  systems  with  comparatively  large  numbers  of  atoms  with  exiting  computers.  An 
example  of  one  such  method  falls  under  the  general  category  of  “atomistic  methods”.  These 
techniques  ASSUME  that  there  is  a  potential  of  interaction  between  atoms.  With  such  an 
assumption,  it  is  possible  to  calculate  dynamic  effects  for  large  ensembles  of  atoms,  e.g. 
modeling  crack-tips,  dislocation  dynamics,  etc.  The  principal  drawback  for  these  methods  is 
that  the  assumed  atomic  potential  can  only  approximate  atomic  interactions. 

When  a  one  assumes  that  this  potential  arises  only  between  pairs  of  atoms,  one  has  a  so  called 
pair  potential.  Pair  potential  techniques  assume  that  the  energy  of  interaction  between  atoms 
depends  only  on  the  distance  between  the  atoms.  There  are  numerous  corrections  to  the 
simple  “pair  approximation”  These  correction  involve  adding  three  and  four  body  terms.  The 
pair  approximation,  by  itself,  is  equivalent  to  assuming  that  only  the  length  or  stretch  of 
chemical  bonds  contributes  to  the  energy  of  the  system.  Addition  of  three  body  terms 
introduces  the  energy  of  bond  bending  to  the  total  energy  of  the  system,  and  higher  order 
terms  introduce  quadrapole  effects  etc..  Because  these  techniques  are  approximate,  it  is  only 
possible  to  calibrate  their  accuracy  against  known  results.  As  an  example,  pair  potential 
techniques  provide  answers  which  are  fairly  good  in  systems  where  there  is  little  directional 
character  to  the  bonding.  So  in  pure  Ni,  Cu,  and  A1  the  calculations  are  reasonable.  However, 
in  Si,  W,  and  Ir  these  calculations  give  poor  results.  Thus  one  is  confronted  with  the 
following  conundrum  when  using  approximate  methods:  while  calculations  work  quite  well  to 
explain  the  properties  of  pure  A1  and  Ni,  what  about  NiAl?  The  only  way  to  tell  is  to  perform 
the  experiment  modeled  and  compare  the  results.  But  if  every  calculation  requires  an 
experiment,  there  is  no  advantage  to  the  calculations. 

Though  some  disagree,  our  opinion  is  that  many  of  these  techniques  are  best  used  to  study 
materials  properties  in  general  rather  than  to  seek  to  model  specific  systems.  A  famous 
example  of  this  kind  of  research  comes  from  bubble  raft  models.  Certainly  these  types  of 
models  provide  much  insight  into  the  dynamics  of  dislocations  in  general,  but  one  would  be 
hard  pressed  to  defend  the  position  that  such  models  have  any  thing  to  do  with  dislocation 
dynamics  in  specific  systems,  e.g.  aluminum. 

There  are  a  number  of  electronic  structure  methods  which  should  be  included  as  approximate 
techniques.  The  first  of  these  is  a  class  of  electronic  structure  methods  known  as  pseudo¬ 
potential  methods.  These  techniques  have  proved  to  be  very  fast,  accurate  and  capable  of 
treating  large  unit  cells  for  certain  types  of  materials.  This  approach  treats  only  the  valence 
electrons  exactly  and  assumes  that  they  see  a  “pseudo-potential”  from  the  core  electrons.  The 
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creation  of  these  pseudo-potentials  is  a  great  art.  Once  a  good  pseudo-potential  is  constructed, 
accurate  calculations  can  be  performed.  The  materials  where  pseudo-potential  techniques 
work  best  are  metals  like  A1  and  Na,  as  well  as  semiconductors  such  as  Si.  A  second 
approximate  method  developed  to  treat  a  dilute  impurity  is  called  the  CPA  or  coherent 
potential  approximation.  Using  this  method  one  does  not  need  to  construct  a  huge  super  cell 
of  100  atoms  to  model  a  1  atomic  per  cent  impurity.  The  CPA  is  a  mean  field  approach  which 
assumes  that  an  atom  in  a  solution  or  alloy  sees  the  average  potential  of  all  the  other  atoms. 
For  many  problems  this  seems  to  be  a  good  approximation;  for  others  it  does  not.  So,  just  as 
with  the  pair  potential  techniques,  because  these  methods  are  approximate  there  can  be  no 
certainty  that  the  results  of  a  calculation  are  predictive  or  simply  manifestations  of  the 
approximation.  In  the  end,  all  approximate  methods  suffer  from  this  uncertainty. 

From  a  metallurgical  perspective  the  biggest  barrier  to  the  use  of  QM  methods  as  practical 
tools  is  the  need  to  treat  larger  cells  and  lower  symmetry  systems.  One  should  keep  in  mind 
however,  that  ten  years  ago  only  those  facilities  with  super-computers  could  undertake  to  do 
QM  calculations  on  crystals  having  more  than  4  atoms  per  cell.  Today,  software  designed  to 
run  on  RISC  workstations  can  be  used  to  calculate  the  electron  density  and  energy  of  crystals 
with  8  or  more  atoms  per  cell.  With  this  pace  of  advancement,  it  is  only  a  matter  of  time 
before  QM  methods  become  traly  useful  as  a  metallurgical  tool.  In  the  interim,  the 
improvement  and  development  of  these  methods  will  be  driven  by  metallurgists  who  use  these 
tools  to  their  present  limit,  and  demand  more  flexible  and  useful  codes.  As  many  of  the 
techniques  are  now  quite  robust,  their  use  by  the  novice  is  constrained  by  the  tedium  of 
creating  input  files  and  analyzing  the  results.  “Builder”  tools  are  now  available  to  set  up  a 
problem  for  analysis  once  the  crystal  structure,  atomic  positions  and  lattice  parameter(s)  are 
known.  Recently  a  number  of  graphic  user  interfaces,  sitting  between  the  user  and  many  of 
these  codes,  have  found  their  way  into  the  marketplace.  What  these  interfaces  lack  is  some 
direction  for  setting  up  metallurgically  meaningful  calculations.  For  example,  the  tools  are 
presently  in  place  which  would  allow  one  to  explore  the  effects  of  ternary  additions  to  the 
phase,  mechanical  behavior,  and  interfacial  energies  of  an  ordered  binary  alloy.  The  design  of 
this  series  of  calculations  would  still  require  an  expert  in  the  use  and  interpretation  of 
electronic  structure  codes.  However,  creating  the  capability,  within  an  existing  GUI,  which 
guides  the  user  through  the  design  of  a  series  of  calculations  and  helps  in  the  analysis  of  the 
results  is  a  feasible  objective.  Such  a  capability  will  only  be  driven  by  demand  and  without 
the  involvement  of  metallurgist  at  this  time,  with  the  tools  that  exist,  this  demand  will  never 
be  realized. 
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3.  EXPERIMENTAL  PROCEDURES 


3.1  Alloy  Casting  and  Homogenization  Heat  Treatment 

Intermetallic  alloys  for  this  program  were  prepared  by  arc-melting.  Arc-casting  was 
accomplished  using  a  specially  constructed,  tri-electrode  apparatus  which  was  designed  to 
accommodate  batch  sizes  up  to  150  grams.  Purified  argon  gas,  gettered  by  flowing  standard 
bottled  gas  over  heated  titanium  chips,  was  used  to  create  a  protective  atmosphere  which 
minimizes  interstitial  pick-up.  As  a  measure  to  reduce  nitrogen  and  oxygen  impurity  levels 
further,  a  titanium  bar  was  melted  in  the  chamber  prior  to  melting  the  alloy  constituents. 
Typical  oxygen  impurity  levels  in  the  argon  gas,  monitored  continuously  throughout  the 
melting  process,  were  below  10’’  ppm  by  weight  (one  part  impurity  in  10*^  parts  argon). 

Buttons  were  cast  in  water-cooled  copper  hearths  using  a  minimum  of  three  melting  and 
resolidification  sequences;  to  minimize  macroscopic  homogeneity,  the  cast  buttons  were 
"flipped"  between  melting  operations. 

Heat  treatments  were  performed  in  an  all-metal  tungsten-element  vacuum  furnace.  Purified 
helium  or  argon  gas,  also  gettered  by  flowing  standard  bottled  gas  over  heated  titanium  chips, 
was  used  to  create  the  protective  atmosphere.  Typical  oxygen  impurity  levels  in  the  inert  gas, 
monitored  continuously  throughout  the  heat  treatment,  were  below  10’’  ppm  by  weight  (one 
part  impurity  in  10*^  parts  inert  gas).  Temperature  was  monitored  by  a  thermocouple  located 
near  the  samples  within  the  hot  zone. 

3.2  Mechanical  Test  Specimen  Preparation 

Specimens  prepared  for  mechanical  testing  were  fabricated  utilizing  a  powder  metallurgy 
process.  The  preparation  of  most  intermetallic  compounds  is  difficult  and  is  attributed  to 
either  volatilization  during  melting,  cracking  during  cool-down  after  solidification,  reaction 
with  apparatus  components  during  hot  pressing  or  HIPing,  or  oxidation  during  heat 
treatments.  These  difficulties  were  considered  before  arriving  at  a  suitable  technique  for 
sample  fabrication. 

The  procedure  selected  for  specimen  preparation  consists  of  (1)  arc-melting  of  materials  to  the 
desired  composition  (2)  crushing,  attriting  and  screening  to  fine  powders  (-325  mesh  or  ^4}x) 
(3)  consolidation  by  vacuum  hot  pressing  and  (4)  EDM  machining  and  mechanically  polishing 
the  surfaces  to  a  600  grit  (<  27|j,)  surface  finish.  Comminution  of  the  arc-melted  buttons  was 
accomplished  using  a  standard  ball  mill  with  tungsten  carbide  balls  under  an  inert  (argon) 
environment. 

This  fabrication  approach  was  used  to  achieve  the  cleanest  possible  material  (low  interstitials 
and  contaminants)  with  the  best  possible  homogeneity.  High  sample  densities  are  required 
during  the  consolidation  step  to  ensure  representative  mechanical  properties.  To  accomplish 
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this,  all  of  the  compositions  were  prepared  by  vacuum  hot  pressing  at  1500-1650°C/70-105 
MPa/1-2  hours.  3.2  cm  X  5.7  cm  graphite  dies  with  molybdenum  foil  liners  and  an  alumina 
release  agent  were  used  for  the  vacuum  hot  pressing. 

3.3  Characterization 

Optical  Microscopy  (OM),  electron  probe  microanalysis  (EPMA),  X-Ray  Diffraction  (XRD), 
and  wet-chemistry  by  inductively  coupled  plasma/atomic  emission  spectroscopy  (ICP)  were 
used  to  characterize  the  structure  and  composition  of  MoSi2  compositions. 

The  microstructure  of  MoSi2  with  the  Cl  Ijj  crystal  structure  was  determined  by  optical  means 
under  polarized  light.  Grain  structures  are  clearly  evident  in  this  imaging  mode.  MoSi2 
samples  with  the  C40  crystal  structure  were  chemically  etched  with  Murakami's  Reagent  (lOg 
K3  Fe[CN]6  +  lOg  KOH  or  NaOH  +  100  ml  H2O)  to  reveal  grain  stmcture. 

Microstructural  phase  analysis  of  various  MoSi2  compositions  were  determined  by  a  Cameca 
MBX  scanning  electron  microprobe  (EPMA)  equipped  with  three  wavelength  spectrometers 
and  a  Kevex  Sigma  energy  dispersive  spectrometer.  Phase  contrast  was  provided  by 
backscattered  electron  imaging.  Compositions  were  determined  by  wavelength  dispersive 
spectroscopy  (WDS)  using  pure  elemental  standards  and  ZAP  corrections. 

Crystal  structures  were  determined  by  x-ray  diffraction  (XRD)  using  a  Rigaku  rotating  anode 
x-ray  system  and  complementary  JCPDS  data  bank  for  crystallographic  identification.  XRD 
with  20  scans  from  20  to  80'^  were  used. 

3.4  Mechanical  Testing 

Specimens  prepared  for  four-point  flexure  testing  exhibited  dimensions  of  approximately  4 
mm  (width)  x  1  mm  (thickness)  x  30  mm  (length).  Testing  was  performed  under  argon  in  an 
all  metal  tungsten-element  vacuum  furnace.  The  molybdenum  fixture  contains  alumina 
contact  rods  with  a  geometry  of  major  span  =  25.4  mm;  minor  span  =  12.7  mm.  Displacement 
was  measured  using  a  center  point  contact  extensometer.  Measurements  of  stress,  modulus, 
and  ductility  were  determined  from  outputted  load-displacement  curves. 

Specimens  for  three-point  Kic  flexure  testing  had  the  same  dimensions  as  the  four-point 
flexure  specimens  with  the  addition  of  a  notch  which  was  mechanically  cut  into  the  sample 
using  a  low  speed  diamond  saw.  The  notch  depth  was  one-half  of  the  specimen  thickness. 


Compression  samples  for  the  determination  of  strength  and  minimum  creep  rates  had 
dimensions  of  approximately  4  mm  x  4  mm  x  10  mm  (height).  Testing  was  performed  in  the 
equipment  described  for  flexural  testing  except  that  SiC  compression  surfaces  were  used. 
Creep  rates  were  determined  from  time-displacement  plots  with  test  durations  of  1(X)0 
seconds  at  each  discreet  increment  of  stress.  One  sample  for  each  temperature  was  used  to 
measure  the  minimum  creep  rate  at  several  stress  levels. 
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4.  RESULTS  AND  DISCUSSION 


4.1  Phase  I  Alloy  Screening 

Six  elemental  modifications  were  selected  for  evaluation.  Aluminum,  boron  and  germanium 
substitutions  for  silicon  and  hafnium,  niobium  and  rhenium  substitutions  for  molybdenum 
were  selected  for  alloy  screening.  These  modifications  will  serve  to  evaluate  the  effects  of 
crystal  structure,  substitutional  alloying,  atomic  radii,  and  e/a  ratio  as  well  as  clarification  or 
confirmation  of  phase  boundaries.  The  initial  samples  for  evaluation  contained  substitutions 
of  50  percent  (16.7  atomic  percent  when  substituted  for  molybdenum  and  33.3  atomic  percent 
when  substituted  for  silicon)  and  were  prepared  by  arc-melting  and  homogenization  heat 
treatment  at  1400°C  for  48  hours.  Additional  samples  of  other  compositions  to  better  define 
phase  boundaries  were  also  prepared  by  arc-melting.  Phase  analyses  were  performed  by  x-ray 
diffraction  (XRD)  and  electron  probe  microanalysis  (EPMA).  The  results  are  discussed 
below. 

4.1.1  Aluminum  Modifications 

The  composition  of  Mo(Sio.5Alo.5)2  exhibited  the  C40  hexagonal  crystal  structure  and  the 
lattice  parameter  is  consistent  with  that  found  by  Nowotny  in  x-ray  card  10-277.  The  XRD 
results  are  summarized  in  Table  3.  Additional  compositions  containing  5  and  40.5  atomic 
percent  aluminum,  Mo(Sio.92Alo.o8)2  and  Mo(Sio.39Alo.6i)2  respectively,  were  also  prepared  to 
better  define  solubilities  and  phase  boundaries.  EPMA  results  generally  confirm  the  phase 
diagram  and  solubility  limits  observed  by  previous  investigators^’.  Backscatter  electron 
images  (BSE)  along  with  the  EPMA  phase  analyses  are  shown  in  Figs.  6,  7  and  8.  The  phase 
boundaries  bounded  by  the  C40  phase  were  found  to  occur  at  aluminum  concentrations  of  8.6 
and  41.6  atomic  percent,  compositions  of  Mo(Sio.87Alo.i3)2  and  Mo(Sio.38Alo.62)2  respectively, 
at  1400°C.  The  solubility  of  aluminum  in  MoSi2  is  3.1  atomic  percent.  Based  on  these 
results,  two  compositions  (20  and  50%  aluminum  substitutions  for  silicon)  of  Mo(Sio.8Alo.2)2 
and  Mo(Sio.5Alo.5)2  were  selected  for  evaluation  of  properties.  Both  of  these  compositions 
have  the  C40  crystal  structure. 

4.1.2  Boron  Modifications 

The  composition  containing  50  percent  substitution  of  boron  for  silicon  [33.3  atomic  percent 
B  or  Mo(Sio.5Bo.5)2]  exhibited  a  multiphase  stmcture  consisting  of  MoSi2  (with  the  Cllb 
stmcture)  and  borides  approximating  the  compositions  MoB  and  M0B2.  A  BSE  image  along 
with  the  EPMA  phase  analysis  is  shown  in  Fig.  9.  These  results  confirm  previous  studies^^ 
and  indicate  no  appreciable  solubility  of  B  in  MoSi2  and  Si  in  "M0B2".  Based  on  these 
results,  boron  modifications  were  not  selected  for  further  evaluation. 
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4. 1 .3  Germanium  Modifications 


The  compositions  containing  10,  25,  and  50  percent  germanium  substitutions  for  silicon  [6.7, 
16.7,  and  33.3  atomic  percent  Ge  or  respective  compositions  of  Mo(Sio.9Geo.i)2, 
Mo(Sio.75Geo.25)2,  and  Mo(Sio.5Geo.5)2]  exhibited  the  single  phase  Cll^  crystal  structure 
indicating  complete  solubility  of  germanium  in  MoSi2  in  the  region  of  interest.  A  BSE  image 
along  with  the  EPMA  phase  analysis  is  shown  in  Fig.  10.  As  shown  by  the  results  in  Table  3, 
the  lattice  parameters  generally  followed  the  rule  of  mixtures  and  increased  with  increasing 
germanium  concentrations.  These  modifications  have  been  selected  for  properties  evaluation. 

4.1.4  Hafnium  Modifications 

The  composition  containing  50  percent  hafnium  substitution  for  molybdenum  [16.7  atomic 
percent  Hf  or  (Moo.5Hfo.5)Si2]  in  MoSi2  exhibited  a  two  phase  structure  of  MoSi2  (with  the 
Cl  Ijj  structure)  and  HfSi2  (with  the  C49  structure).  This  result  verifies  phase  diagram  results 
by  other  investigators^^  and  indicates  very  low  solubilities  of  Hf  in  MoSi2  and  Mo  in  HfSi2.  A 
BSE  image  along  with  the  EPMA  phase  analysis  is  shown  in  Fig.  11.  Based  on  these  results, 
hafnium  modifications  were  not  selected  for  further  evaluation. 

4.1.5  Niobium  Modifications 

The  composition  containing  a  50  percent  niobium  substitution  for  molybdenum  [16.7  atomic 
percent  Nb  or  (Moo.5Nbo.5)Si2]  exhibited  the  C40  hexagonal  crystal  structure.  An  additional 
composition  containing  6.3  atomic  percent  Nb  [(Moo.8iNbo.i9)Si2]was  prepared  to  better 
define  the  solubility  limits.  The  results  of  EPMA  analyses  indicate  solubility  limits  that  are 
different  than  that  observed  by  other  investigators^''.  BSE  images  along  with  the  EPMA  phase 
analyses  are  shown  in  Figs.  12  and  13.  The  solubility  of  niobium  in  MoS^  (Cl  Ij,  structure)  is 

1.5  atomic  percent,  (Moo.95Nbo.o5)Si2, while  the  solubility  of  molybdenum  in  NbSi2  (C40 
structure)  is  26.5  atomic  percent,  (Moo.8Nbo.2)Si2 ,  at  1400°C.  XRD  results  in  Table  3  indicate 
the  lattice  parameter  of  MoSi2  (C40  structure)  increases  linearly  with  increasing  niobium 
concentrations.  Based  on  these  results,  compositions  with  9.7  and  16.7  atomic  percent 
niobium  substitutions  for  molybdenum,  (Moo.79Nbo.29)Si2  and  (Moo.5Nbo,5)Si2,  in  MoSi2  were 
selected  for  further  evaluation.  Both  of  these  alloys  lie  within  the  C40  phase  field. 

4.1.6  Rhenium  Modifications 

The  compositions  containing  10,  25,  and  50  percent  rhenium  substitutions  for  molybdenum 
[3.3,  8.3,  and  16.7  atomic  percent  Re  or  (Moo.9Reo.i)Si2,  (Moo.75Reo.25)Si2,  and  (Moo.5Reo.5)Si2 
respectively]  exhibited  the  Cl  1^  tetragonal  structure.  A  typical  BSE  image  along  with  the 
EPMA  phase  analysis  is  shown  in  Fig.  14.  The  results  indicate  complete  solubility  in  the 
ReSi2-MoSi2  system.  XRD  results  in  Table  3  indicate  the  lattice  parameter  decreases  linearly 
with  increasing  rhenium  concentrations.  These  compositions  were  selected  for  further 
evaluation. 
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4.2  Phase  I  Property  Evaluation 

4.2.1  Mechanical  Properties 

The  mechanical  properties  have  been  determined  by  four  point  flexural  testing.  The  testing 
was  performed  at  room  temperature  and  at  elevated  temperatures  to  determine  strength, 
ductility,  and  ductile-to-brittle  transition  temperature  (DBTT).  These  results  are  summarized 
in  Table  4.  In  general,  no  ductility  was  observed  below  the  DBTT.  The  strength  of  the  alloys 
appear  to  increase  with  increasing  temperatures  while  the  elastic  modulus  decreases  with 
increasing  temperature.  The  DBTT  was  not  altered  significantly  with  the  alloy  compositions 
and  crystal  structures  evaluated.  The  observed  DBTT  appears  to  be  higher  than  that  reported 
in  the  literature"*  and  may  be  attributed  to  the  processing  techniques  used.  These  samples 
exhibit  high  purities  and  low  interstitial  contents  which  may  explain  the  observed  higher 
DBTTs. 

4.2.2  Oxidation  Resistance 

Samples  for  the  evaluation  of  oxidation  resistance  and  mechanical  properties  were  prepared  in 
the  same  manner.  A  procedure  of  arc-melting/powder  comminutionA^HP  consolidation  was 
utilized.  Oxidation  testing  was  performed  at  1400°C  for  22  hours  and  124  hours  under 
isothermal  conditions.  This  test  was  a  screening  test  to  determine  if  the  alloy  modifications 
selected  would  degrade  oxidation  resistance. 

The  weight  change  results  shown  in  Table  5  indicate  none  of  the  alloy  modifications 
significantly  decreased  the  oxidation  resistance  of  MoSia.  However,  the  samples  containing 
50%  aluminum  substitutions  exhibited  spalling  (apparently  due  to  the  formation  of  a  mixed 
oxide  scale)  while  the  composition  containing  50%  substitution  of  rhenium  for  molybdenum 
exhibited  large  initial  weight  losses. 

4.3  Phase  II  Alloy  Screening 

Six  alloying  elements  were  evaluated  in  the  Phase  I  effort.  The  ternary  alloy  modifications  of 
MoSia  include  Al,  B,  Ge,  Hf,  Nb,  and  Re.  Based  on  the  results  in  Phase  I  which  concluded 
that  B  and  Hf  have  virtually  no  solubility  in  MoSia,  these  modifications  were  not  selected  for 
further  evaluation. 

Four  alloying  elements  were  selected  for  the  evaluation  of  mechanical  properties:  Aluminum 
and  germanium  substitutions  for  silicon  in  MoSia  and  niobium  and  rhenium  substitutions  for 
molybdenum  in  MoSia.  All  alloy  compositions  selected  in  this  phase  of  the  program  are 
single  phase.  Alloy  compositions  for  the  last  phase  of  this  program  were  either  near  the  Cl  1^ 
/C40  two  phase  boundaries  or  within  the  two  phase  field.  These  modifications  were  served  to 
evaluate  the  effects  of  crystal  structure,  substitutional  alloying,  atomic  radii  size,  and  e/a  ratio 
on  mechanical  properties.  Samples  of  the  selected  compositions  were  prepared  by  arc- 
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melting,  comminution,  blending,  and  vacuum  hot  pressing.  Samples  prepared  using  this 
methodology  have  microstructures  typified  by  the  baseline  MoSia  microstructure  shown  in 
Fig.  15. 

4.3.1  Aluminum  Modifications 

The  compositions  of  Mo(Sio.8Alo.2)2  and  Mo(Sio.5Alo.5)2  exhibiting  the  C40  hexagonal  crystal 
structure  were  selected  for  properties  evaluation.  These  compositions  both  lie  within  the  C40 
phase  boundaries  defined  previously  at  Mo(Sio.87Alo.i3)2  and  M0AI2. 

4.3.2  Germanium  Modifications 

The  compositions  Mo(Sio.9Geo.i)2  and  Mo(Sio.75Geo.25)2  exhibiting  a  single  phase  Cllb 
crystal  structure  with  complete  solubility  of  germanium  in  MoSia  were  selected  for  properties 
evaluation. 

4.3.3  Niobium  Modifications 

The  compositions  (Moo.7iNbo.29)Si2  and  (Moo.5Nbo.5)Si2  exhibiting  the  C40  hexagonal  crystal 
structure  were  selected  for  properties  evaluation.  These  compositions  both  lie  within  the  C40 
phase  boundaries  defined  previously  at  (M00.8  Nbo.2)Si2.and  NbSia. 

4.3.4  Rhenium  Modifications 

The  compositions  (Moo^Reo.OSia,  (Moo.75Reo.25)Si2,  (Moo.5Reo.5)Si2  exhibiting  the  Cllb 
tetragonal  structure  with  complete  solubility  of  rhenium  in  MoSia  were  selected  for  properties 
evaluation. 

4.4  Phase  II  Properties  Evaluation 

4.4.1  Mechanical  Properties 

4.4. 1 . 1  Four-Point  Flexural  Testing 

The  mechanical  properties  for  the  baseline  MoSia  and  modifications  with  Al,  Ge,  Nb,  and  Re 
have  been  determined  by  four  point  flexural  testing.  The  testing  was  performed  at  room 
temperature  and  at  elevated  temperatures  to  determine  strength,  ductility,  and  ductile-to-brittle 
transition  temperature  (DBTT).  These  results  are  summarized  in  Table  6.  In  general,  no 
ductility  was  observed  below  the  DBTT.  The  strength  of  the  alloys  appear  to  increase  with 
increasing  temperatures  while  the  elastic  modulus  decreases  with  increasing  temperature.  The 
DBTT  was  not  altered  significantly  with  the  alloy  compositions  and  crystal  structures 
evaluated.  The  observed  DBTT  were  all  in  the  range  1250  to  1350°C. 
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The  results  in  Table  6  also  indicate  that  strengthening  was  significant  with  the  niobium  alloy 
modification  at  elevated  temperatures. 

4.4. 1 .2  Three-Point  Flexural  Kjc  Testing 

The  fracture  toughness  was  determined  using  a  notched  specimen  tested  in  three  point  flexure. 
The  results  for  the  baseline  MoSi2  and  the  Al,  Ge,  Nb,  and  Re  modified  alloys  are  shown  in 
Table  7.  Generally,  all  of  the  alloys  have  low  fracture  toughnesses  below  the  DBTT  (1-4  MPa 
•Vm)  and  improved  fracture  toughnesses  when  tested  above  the  DBTT. 

The  niobium  modified  alloy  exhibited  a  higher  fracture  toughness  at  elevated  temperatures 
than  either  the  baseline  MoSi2  or  the  aluminum,  germanium,  or  rhenium  modified  alloys. 
This  is  most  likely  attributed  to  its  higher  elevated  temperature  strength. 

4.4. 1.3  Compression  Testing 

Compression  testing  was  performed  to  determine  suitable  stresses  to  perform  compressive 
creep  testing.  A  sample  of  each  composition  was  used  to  determine  the  0.1%  offset  yield 
stress  throughout  the  temperature  range  1000-1400°C.  The  first  test  was  performed  at  1400°C 
to  a  strain  of  0.1%.  Subsequent  tests  were  performed  on  the  same  sample  to  a  similar  strain 
level  with  decrements  of  100°C. 

These  results  are  shown  in  Table  8  and  Fig.  16  and  indicate  that  the  addition  of  niobium 
dramatically  improves  the  strength  of  MoSia  throughout  the  temperature  range.  Significant 
strengthening  is  also  observed  with  aluminum  modifications  at  the  lower  temperatures. 
Germanium  modifications  were  found  to  reduce  the  strength  of  MoSia.  Rhenium 
modifications  appear  to  increase  the  strength  of  MoSi2  based  on  the  limited  data  generated. 
Difficulties  encountered  during  testing  limited  the  amount  of  fully  reduced  data  for  the 
rhenium  modified  alloys. 

4.4. 1.4  Compressive  Creep  Testing 

Compressive  creep  testing  was  performed  for  the  baseline  MoSi2  and  the  Al,  Ge,  Nb,  and  Re 
modifications  to  determine  the  minimum  creep  rate  at  temperatures  of  1000,  1200,  and  1400° 
C.  One  test  specimen  was  used  at  each  temperature  and  creep  rates  were  determined  starting 
at  the  lowest  test  stresses.  Loading  was  maintained  at  each  stress  level  for  1000  seconds  to 
determine  minimum  creep  rates.  An  approximation  of  the  creep  ductility  was  determined  by 
totaling  the  creep  strains  for  all  of  the  applied  stress  levels  at  one  test  temperature. 

The  minimum  creep  rates  for  each  composition  are  listed  in  Tables  9, 10,  and  1 1  and  shown  in 
Figs.  17,  18,  and  19  for  test  temperatures  of  1000,  1200,  and  1400°C  respectively.  The  total 
creep  strain  at  each  test  temperature  is  tabulated  in  Table  12.  These  results  indicate  that 
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modifications  with  germanium  generally  decrease  creep  resistance.  The  addition  of  aluminum 
improved  creep  strength  only  at  the  lower  temperatures.  The  addition  of  rhenium  improved 
creep  resistance  by  about  an  order  of  magnitude  and  is  consistent  with  prior  work^.  However, 
this  improvement  in  creep  resistance  is  accompanied  by  a  significant  reduction  in  the  creep 
ductility  of  the  alloy.  The  addition  of  niobium  to  MoSi2  also  provided  a  significant 
improvement  in  the  creep  resistance  by  about  an  order  of  magnitude  throughout  the 
temperature  range  of  interest  while  not  debiting  the  creep  ductility  of  the  alloy.  The  need  to 
improve  the  high  temperature  strength  of  MoSii  appears  to  have  been  met  with  the  additions 
of  Re  and  Nb.  The  addition  of  Nb  offers  the  best  promise  since  compressive  strength  and 
creep  strength  were  increased  significantly  without  a  debit  in  alloy  ductility. 

The  stress  exponent  was  calculated  from  the  data  in  Figs.  17,  18,  and  19.  Generally,  n=3-5  at 
1000°C,  n=3  at  1200°C,  and  n=2  at  1400°C.  These  n  values  are  in  the  “power  law  creep” 
region  where  vacancy  diffusion  assisted  dislocation  climb  occurs. 

4.5  Phase  III  Alloy  Screening 

Six  alloying  elements  were  in  Phase  I  of  this  program.  The  ternary  alloy  modifications  of 
MoSi2  include  Al,  B,  Ge,  Hf,  Nb,  and  Re.  Based  on  this  work,  Al,  Ge,  Nb  and  Re 
modifications  were  selected  for  further  evaluation  of  properties. 

Four  alloying  elements  were  selected  for  evaluation  of  mechanical  properties  in  the  Phase  n 
effort:  Aluminum  and  germanium  substitutions  for  silicon  in  MoSi2  and  niobium  and 
rhenium  substitutions  for  molybdenum  in  MoSi2  were  selected.  All  alloy  compositions 
selected  in  this  phase  of  the  program  were  single  phase.  The  mechanical  properties  testing 
during  the  second  year  effort  indicated  improvements  in  creep  strength  imparted  by  additions 
of  either  niobium  or  rhenium^^. 

Alloy  compositions  for  this  phase  of  the  program  are  niobium  and/or  rhenium  modifications 
near  the  Cl  1|3/C40  two  phase  boundaries  or  within  the  two  phase  field.  These  modifications 
serve  to  evaluate  the  effects  of  crystal  structure,  substitutional  alloying,  atomic  radii  size,  and 
e/a  ratio  on  mechanical  properties.  Samples  of  the  selected  compositions  were  prepared  by 
the  standard  arc-melting,  comminution,  blending,  and  vacuum  hot  pressing  used  in  Phases  I 
and  n. 

4.5.1  Niobium  Modifications 

The  compositions  (M00.97  Nbo.o3)Si2  with  the  Cl  lb  tetragonal  crystal  structure;  (M00.93 

Nbo.o7)Si2  and  (Moo.gs  Nbo.i5)Si2  with  the  two  phase  C40  +  Cl  lb  structure;  and  (M00.71 
Nbo.29)Si2  exhibiting  the  C40  hexagonal  crystal  structure  were  selected  for  properties 
evaluation. 
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4.5.2  Rhenium  Modifications 


The  composition  (Moo.9Reo.i)Si2  exhibiting  the  Cllj)  tetragonal  structure  with  complete 
solubility  of  rhenium  in  MoSi2  was  selected  for  properties  evaluation.  The  mechanical 
properties  measured  during  the  second  year  indicated  the  low  rhenium  modification  exhibited 
the  best  combination  of  mechanical  properties  and  oxidation  resistance. 

4.5.3  Niobium  and  Rhenium 

An  alloy  incorporating  both  niobium  and  rhenium  was  prepared  with  the  composition 
(Moo,6iNbo.29Reo.io)Si2.  This  single  phase  composition  was  selected  based  on  the  excellent 
properties  achieved  with  the  niobium  and  rhenium  alloy  modifications. 

4.6  Phase  III  Properties  Evaluation 

4.6.1  Mechanical  Properties 

4.6. 1 . 1  Four-Point  Flexural  Testing 

The  mechanical  properties  for  the  Nb,  Re  and  Nb  +  Re  modifications  have  been  determined 
by  four  point  flexural  testing.  The  testing  was  performed  at  room  temperature  and  at  elevated 
temperatures  to  determine  strength,  ductility,  and  ductile-to-brittle  transition  temperature 
(DBTT).  These  results  are  summarized  in  Table  13.  In  general,  no  ductility  was  observed 
below  the  DBTT.  The  strength  of  the  alloys  appear  to  increase  with  increasing  temperatures 
while  the  elastic  modulus  decreases  with  increasing  temperature.  The  DBTT  was  not  altered 
significantly  with  the  alloy  compositions  and  crystal  structures  evaluated.  The  observed 
DBTT  appears  to  be  higher  than  that  reported  in  the  literature'*  and  may  be  attributed  to  the 
processing  techniques  used.  Higher  purities  with  lower  interstitial  containing  samples  were 
evaluated  under  this  program. 

The  results  in  Table  13  also  indicate  that  strengthening  was  significant  with  the  niobium  alloy 
modification  at  elevated  temperatures. 

4.6. 1.2  Three-Point  Flexural  Kjc  Testing 

The  fracture  toughness  was  determined  using  a  notched  specimen  tested  in  three  point  flexure. 
The  results  for  the  Nb,  Re  and  Nb  +  Re  modified  MoSi2  are  shown  in  Table  14.  Generally,  all 
of  the  alloys  have  low  fracture  toughnesses  (1-4  MPa  Vm)  below  the  DBTT  and  improved 
fracture  toughnesses  when  tested  above  the  DBTT. 

The  niobium  and/or  rhenium  modified  alloys  containing  9.5  atomic  percent  Nb  exhibited  a 
slightly  higher  fracture  toughness  at  elevated  temperatures  than  either  the  other  modified 
alloys. 


4.6. 1.3  Compression  Testing 

Compression  testing  was  performed  to  determine  suitable  stresses  to  perform  compressive 
creep  testing.  A  sample  of  each  composition  was  used  to  determine  the  0.1%  offset  yield 
stress  throughout  the  temperature  range  1000-1400°C.  The  first  test  was  performed  at  1400°C 
to  a  strain  of  0.1%.  Subsequent  tests  were  performed  on  the  same  sample  to  a  similar  strain 
level  with  decrements  of  100°C. 

The  results  are  shown  in  Table  15  and  Fig.  20  and  indicate  the  addition  of  niobium,  rhenium, 
and  niobium  +  rhenium  improves  the  strength  of  MoSi2  throughout  the  temperature  range. 

4.6. 1.4  Compressive  Creep  Testing 

Compressive  creep  testing  was  performed  for  the  Nb,  Re,  and  Nb  +  Re  modifications  to 
determine  the  minimum  creep  rate  at  temperatures  of  1000,  1200,  and  1400°C.  One  test 
specimen  was  used  at  each  temperature  and  creep  rates  were  determined  starting  at  the  lowest 
test  stresses.  Loading  was  maintained  at  each  stress  level  for  1000  seconds  to  determine 
minimum  creep  rates.  An  approximation  of  the  creep  ductility  was  determined  by  totaling  the 
creep  strains  for  all  of  the  applied  stress  levels  at  one  test  temperature. 

The  minimum  creep  rates  for  each  composition  are  listed  in  Tables  16,  17,  and  18  and  shown 
in  Figs.  21,  22,  and  23  for  test  temperatures  of  1000,  1200,  and  1400°C  respectively.  The 
total  creep  strain  at  each  test  temperature  is  tabulated  in  Table  19.  The  addition  of  rhenium, 
niobium,  and  Nb  +  Re  improved  creep  resistance  by  an  order  of.  However,  this  improvement 
in  creep  resistance  is  accompanied  by  a  significant  reduction  in  the  creep  ductility  of  the  Re 
modified  alloy.  The  addition  of  niobium  and  Nb  +  Re  to  MoSi2  provided  a  significant 
improvement  in  the  creep  resistance  while  not  debiting  the  creep  ductility  of  the  alloy. 

The  stress  exponent  was  calculated  from  the  data  in  Figs.  21,  22,  and  23.  Generally,  n=3-7  at 
1000°C,  n=3  at  1200°C,  and  n=2  at  1400°C.  These  n  values  are  in  the  “power  law  creep” 
region  where  vacancy  diffusion  assisted  dislocation  climb  occurs  except  at  1000°C.  The 
stress  exponent  n  =  5-7  for  several  alloys  when  tested  at  the  highest  stress  levels  (about  1000 
Mpa)  at  1000°C  suggests  a  mechanism  where  thermally  activated  dislocation  glide  is  the  rate¬ 
controlling  step,  perhaps  associated  with  a  significant  back  stress  or  threshold  stress. 

4.6.2  Oxidation  Testing 

Samples  for  evaluation  of  oxidation  resistance  were  prepared  in  the  same  manner  as  the 
samples  used  for  mechanical  properties  evaluation.  Samples  were  polished  to  a  600  grit  finish 
prior  to  oxidation  testing.  They  were  tested  in  air  at  1400°C  for  20  hours  and  120  hours 
under  isothermal  conditions.  This  test  was  to  determine  if  the  alloy  modifications  selected 
would  degrade  oxidation  resistance. 
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The  weight  change  results  shown  in  Table  20  indicate  all  of  the  alloy  modifications  did  not 
debit  the  excellent  oxidation  resistance  of  MoSia.  The  MoSi2  weight  change  results  are 
included  from  the  first  oxidation  test  for  comparison  purposes. 

4.7  Atomic  Modeling 

A  series  of  atomic  modeling  calculations  was  done  for  MoSi2  to  determine  the  feasibility  of 
using  efficient  simulation  tools  to  assist  in  the  interpretation  of  properties,  specifically 
bonding  and  the  influence  of  alloying.  This  work  was  mostly  carried  out  with  internal  funds, 
and  was  facilitated  by  Mark  Eberhart  of  the  Colorado  School  of  Mines. 

The  LKKR^^  code  had  been  installed  at  UTRC  and  was  successfully  used  to  model  antiphase 
boundaries  in  NisAl.  The  baseline  Cl  lb  structure  is  shown  in  Figure  24.  An  input  file  was 
created  for  MoSi2  and  the  potentials  were  converged  on  an  Alliant  minisupercomputer.  The 
charge  density  was  then  computed.  These  calculations  took  about  one  day.  For  visualization 
purposes,  the  charge  density  was  converted  from  a  2D  to  a  3D  structure  in  the  “PLOT3D” 
format  for  mesh  co-ordinates  and  relevant  quantities.  To  avoid  heavy  contouring  near  the 
nucleus,  another  routine  was  used  to  crop  the  data;  this  allows  one  to  focus  on  the  bonding 
electrons.  It  was  found  that  the  gradient  of  charge  density  was  useful  for  locating  bonds 
within  a  structure.  A  routine  was  written  and  used  to  compute  the  x,  y  and  z  components  of 
the  charge  density. 

The  UTRC  visualizers  “VISA”  and  “VISA-GOLD”  were  used  to  examine  the  charge  density. 
A  related  commercial  visualizer,  “FIELD VIEW”  was  also  used.  An  example  of  a  set  of 
charge  density  contours  on  the  (001)  plane  is  shown  in  Figure  25.  It  can  be  seen  that  the 
expected  four-fold  symmetry  is  not  completely  satisfied.  However,  some  indication  of 
elongation  of  the  charge  density  around  Si  atoms  towards  the  central  Mo  atom  can  be  seen  for 
the  (200)  section  in  Figure  26.  Three  dimensional  viewing  of  the  charge  density  using 
“isosurfaces”  showed  sharp  spikes  in  the  region  where  bonds  were  expected.  It  was 
concluded  that  these  were  associated  with  numerical  noise  and  that  the  LKKR  code  was  not 
adequate  to  carry  out  precise  calculations  at  the  low  charge  densities  in  the  interatomic 
regions. 

The  LASTO  code  was  subsequently  installed  at  UTRC  to  run  on  an  IBM  RISC  6000 
workstation  using  the  “Distributed  Queuing  System”  tool  for  distributed  computing  on  a 
network  of  about  a  dozen  such  workstations.  The  most  available  workstation  was  selected 
and  a  single  processor  was  used.  Scripts  were  built  to  compile,  time  and  track  the  submitted 
jobs.  Machines  with  128MB  of  RAM  were  used.  Data  were  stored  on  a  DECstation  running 
Ultrix,  or  on  an  Epoch  archival  file  server.  A  “Builder”  was  installed  to  facilitate  the  building 
of  molecular  models.  Symmetry  files  and  the  related  input  files  for  structure  and  initial 
potentials  were  built  and  tested  for  A1  and  Ag  (fee);  NiAl  and  CoAl  (B2);  NiaAl  (LI2);  TisAl 
(Llo);  Fe  (bcc);  Co  (hep);  and  MoSi2  (Cl  lb). 
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Again,  the  potentials  were  converged  and  then  charge  density  distributions  were  computed  at 
various  resolutions.  Tools  were  used  to  convert  to  PLOT3D,  crop  and  compute  gradients. 
These  calculations  took  about  one  day  of  CPU  time.  The  FIELDVIEW  isosurface  and 
“cutting  plane”  tools  were  used  to  carry  out  viewing.  A  SGI  indigo2  extreme  workstation  was 
used  to  view  the  charge  density.  The  “snapshot”  and  “showcase”  tools  were  used  to  create, 
size  and  convert  images  into  the  Adobe  PostScript  format  for  printing  as  black  &  white  or 
color  figures.  A  “surface”  (filled)  view  of  the  charge  density  on  the  (110)  plane  is  shown  in 
Figure  27.  The  Mo  nuclei  can  clearly  be  seen  at  the  center  and  the  comers  in  this  plot.  The 
other  six  atom  centers  are  Si. 

These  results  were  generally  much  better  behaved  (smoother)  than  the  LKKR  calculations.  As 
a  result,  contour  plots,  such  as  seen  in  Figures  28  and  29  show  elongation  of  the  Si  electrons 
toward  the  Mo.  The  synunetry  is  very  close  to  what  is  expected.  Bonds  can  be  clearly  seen  in 
Figure  28.  They  are  in  the  expected  directions.  Assessment  of  the  strength  of  the  bonds  is 
subject  to  interpretation  at  this  point,  but  it  can  be  said  that  this  material  has  a  strong  covalent 
component  since  the  electrons  are  concentrated  along  the  bond  lines.  This  would  not  be  the 
case  for  a  metal,  such  as  Ni  or  even  an  intermetallic  that  behaves  like  a  metal,  e.g.  NisAl.  A 
full  cell  view  of  the  charge  density  at  a  value  where  the  contours  overlap  significantly  is 
shown  in  Figure  30.  The  distortion  of  the  charge  around  the  silicon  atoms  toward  the 
molybdenum  sites  is  clearly  evident. 

Such  figures  can  be  used  to  rationalize  cleavage  behavior,  but  guidance  for  alloy  design  is  not 
immediately  obvious.  Attempts  to  add  third  elements  to  such  a  calculation  for  the  case  of 
Ni3Al  are  now  underway;  however,  this  is  challenging  on  a  single  processor  workstation,  at 
least  for  anything  other  than  a  major  addition.  A  more  efficient  version  of  LASTO  is  about  to 
be  released  by  James  McLaren  of  Tulane  University.  This  upgraded  software  and  continual 
advances  in  harware  capability  should  assist  in  carrying  out  such  challenging  calculations. 
Similar  tools  are  routinely  used  in  the  petrochemical,  pharmaceutical  and  semiconductor 
industries,  and  it  is  only  a  matter  of  time  before  they  are  used  by  the  Materials  Community  as 
well. 
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5.  CONCLUSIONS 


1.  Six  elemental  modifications  of  MoSi2  have  been  evaluated.  They  are  Al,  B,  and  Ge 
subsitutions  for  Si  and  Hf,  Nb,  and  Re  subsitutions  for  Mo. 

2.  B  and  Hf  modifications  exhibited  two  phase  or  multiphase  structures  with  very  little 
solubility  in  MoSi2. 

3.  Al  and  Nb  substitutions  for  Si  and  Mo,  respectively,  in  MoSi2  resulted  in  alloys  with  the 
C40  hexagonal  crystal  structure. 

4.  Ge  and  Re  substitutions  for  Si  and  Mo,  respectively,  in  MoSi2  resulted  in  alloys  with  the 
Cl  I5  tetragonal  crystal  structure. 

5.  Most  of  the  analyses  confirm  the  validity  of  available  phase  diagrams  although  the  results 
presented  here  indicate  a  need  to  refine  the  solubility  limits. 

6.  Isothermal  oxidation  testing  indicates  no  significant  debit  in  oxidation  resistance  at  1400° 
C  due  to  ternary  alloying.  Alloys  containing  higher  concentrations  of  ternary  elements 
may  exhibit  reduced  oxidation  resistances. 

7.  Mechanical  properties  as  determined  by  four-point  flexural  testing  indicate  DBTTs  were 
in  the  temperature  range  1250-1 350°C  for  all  of  the  modifications  evaluated. 

8.  The  fracture  toughness  of  the  baseline  and  all  of  the  modified  alloys  were  1-4  MPa  Vm 
below  the  DBTT  and  slightly  higher  (l-13MPa  Vm)  above  the  DBTT.  The  9.5  atomic 
percent  niobium  alloys  with  the  C40  structure  exhibited  a  better  fracture  toughness  (5-13 
MPa-Vm)  at  elevated  temperatures. 

9.  The  compressive  strength  of  MoSi2  was  improved  significantly  with  the  addition  of  Al, 
Nb  and  Re.  Ge  did  not  improve  the  mechanical  properties  of  MoSi2.while  Re  improved 
the  elevated  temperature  strength  The  strength  increase  due  to  Al  was  not  retained  at  the 
highest  test  temperature  (1400°C).  Nb  and  Re  additions  increased  strength  at  all 
temperatures  tested.. 

10.  The  compressive  minimum  creep  rate  was  improved  with  the  additions  of  Nb  and  Re.  The 
dramatic  reduction  in  the  minimum  creep  rate  of  the  Re  modified  alloys  was  accompanied 
by  a  reduction  in  the  creep  ductility.  Nb  modified  alloys  also  exhibited  a  significant 
reduction  in  the  compressive  minimum  creep  rate  and  a  reduction  in  creep  ductility  was 
not  observed  when  compared  to  the  baseline  MoSi2. 
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TABLE  1  ATOMIC  RADII 


Element 

Goldschmidt  Radius  (nm) 

Si 

0.117 

B 

0.097 

C 

0.077 

A1 

0.143 

Ge 

0.139 

Mo 

0.140 

Ti 

0.147 

V 

0.136 

Cr 

0.128 

Zr 

0.160 

Nb 

0.147 

Hf 

0.159 

Ta 

0.147 

W 

0.141 

Re 

0.138 
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TABLE  2  WET  CHEMICAL  ANALYSES 


Weight  Percent 


Mo 

Nb 

Re 

Si 

A1 

Ge 

Nominal 

63.04 

36.96 

Actual 

66.16 

35.49 

Mo(Sio.5Alo.5)2 

Nominal 

63.60 

18.56 

17.83 

Actual 

63.48 

17.52 

17,52 

Mo(Sio.9Geo.i)2 

Nominal 

59.53 

31.40 

9.06 

Actual 

62.75 

30.70 

8.14 

Mo(Sio.7sGeo.25)2 

Nominal 

54.98 

24.16 

20.86 

Actual 

57.17 

23.97 

18.62 

(Moo.5Nbo.5)Si2 

Nominal 

31.89 

30.88 

37.23 

Actual 

31.85 

29.15 

34.87 

Nominal 

53.64 

11.45 

34.91 

Actual 

56.48 

11.10 

33.38 

(MOo.75RCo.2s)Si2 

Nominal 

41.23 

26.57 

32.20 

Actual 

43.19 

25.90 

30.76 

(MOo.5RCo.5)Si2 

Nominal 

24.34 

47.24 

28.42 

Actual 

24.94 

46.77 

26.81 

Lattice  6-681  VHP  VHP 

Plane  MoSi2 _ MoSU _ (Moo.9Reo.i)Si2 


6.558  6.592 


TABLE  4  FOUR-POINT  FLEXURE  TEST  RESULTS 


TABLE  5  ISOTHERMAL  OXIDATION  RESULTS  (1400°C) 


TABLE  6  FOUR-POINT  FLEXURE  TEST  RESULTS 


36 
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TABLE  8  COMPRESSIVE  YIELD  STRESS  TEST  RESULTS 


0.1%  Yield  Stress  (MPa) 


Composition 

1000°C 

1100°C 

1200°C 

1300°C 

1400°C 

799 

661 

455 

365 

165 

1213 

937 

503 

296 

158 

Mo(Sio.9Geo.i)2 

710 

579 

303 

179 

117 

Mo(Sio.75Geo.25)2 

441 

365 

207 

124 

83 

(MOo.7lNbo.29)Si2 

1158 

834 

737 

489 

351 

(Moo.9Reo.i)Si2 

>1013 

>331 

>145 

(Moo.5Reo.5)Si2 

>937 

>372 

186 

TABLE  9  MINIMUM  CREEP  RATES  AT  1000°C 


13.3  A1 

9.5  Nb 

6.7  Ge 

16.7  Ge 

3.3  Re 

16.7  Re 

TABLE  10  MINIMUM  CREEP  RATES  AT  1200°C 


Stress  (MPa) 


14 

28 

41 

55 

69 

83 

97 

110 

124 

152 

179 

207 

221 

RATES  AT  1400°C 


[luin  Creep  Rate  (s'^) 


6.7  Ge 

16.7  Ge 

3.3  Re 

16.7  Re 

1.7E-06 

3.4E-07 

l.OE-06 

2.0E-06 

5.0E-06 

2.8E-06 

1.2E-06 

6.0E-06 

l.lE-05 

1.6E-05 

l.lE-05 

2.5E-05 

2.7E-06 

2.5E-05 

3.3E-06 

4.5E-06 

6.5E-06 

1 

TABLE  12  COMPRESSIVE  CREEP  DUCTILITY  RESULTS 


Creep  Ductility  (%) 


Composition 

Creep  Region  of  interest 

1000°C 

1200°C 

1400°C 

i+n 

2.0 

4.7 

18.9 

i+n+m 

8.5 

16.4 

29.1 

Mo(Sio.8AIoa)2 

i+n 

0.6 

3.2 

7.5 

i+n+m 

1.4 

9.3 

18.4 

^bhihhi 

I+n 

0.6 

2.9 

2.3 

I+n+m 

2.6 

13.9 

19.4 

Mo(Sio.75Geo.25)2 

I+n 

0.6 

0.3 

0.5 

I+n+m 

2.9 

8.3 

8.3 

I+n 

1.0 

2.8 

7.5 

I+n+m 

1.5 

9.0 

15.5 

(Moo.9Reo.i)Si2 

I+n 

0.6 

0.8 

2.7 

I+n+m 

1.0 

2.6 

6.0 

(MOo.5RC0.5)Si2 

I+n 

0.7 

3.3 

1.9 

I+n+m 

0.7 

5.2 

9.3 
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TABLE  13  FOUR-POINT  FLEXURE  TEST  RESULTS 


TABLE  14  THREE-POINT  FLEXURE  Kic  TEST 
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TABLE  15  COMPRESSIVE  YIELD  STRESS  TEST  RESULTS 


0.1%  Yield  Stress  (MPa 


Composition 

1000°C 

1300°C 

1400°C 

799 

661 

455 

365 

165 

(MOo.97Nbo.03)Si2 

924 

731 

469 

309 

162 

(MOo.93  Nbo.07)Si2 

814 

658 

485 

315 

181 

(MOo.85  Nbo.l5)Si2 

1186 

917 

654 

450 

283 

(MOo.71  Nbo.29)Si2 

1062 

786 

519 

319 

188 

(Moo.9Reo.i)Si2 

1234 

,  758 

450 

318 

181 

(Moo.6iNbo.29Reo.i)Si2 

1579 

1065 

580 

261 

134 

■6 


TABLE  16  MINIMUM  CREEP  RATES  AT  1000°C 


Minimum  Creep  Rate  (s‘^) 


Stress  (MPa) 

INb 

2.5  Nb 

5Nb 

9.5  Nb 

3.3  Re 

9.5Nb-3Re 

84.5 

8.0E-7 

4.2E-7 

85.7 

6.1E-7 

2.8E-7 

160 

1.9E-6 

169 

1.9E-6 

1.4E-6 

171 

1.4E-6 

8.6E-7 

l.lE-6 

240 

3.8E-6 

254 

4.4E-6 

4.1E-6 

2.6E-6 

257 

4.8E-6 

1.3E-6 

321 

7.9E-6 

338 

8.9E-6 

1.9E-5 

343 

1.3E-5 

1.5E-6 

4.5E-6 

401 

1.7E-5 

423 

1.8E-5 

429 

3.8E-5 

3.0E-6 

8.0E-6 

481 

4.3E-5 

507 

3.5E-5 

514 

1.2E-4 

4.9E-6 

1.5E-5 

592 

7.0E-5 

600 

8.3E-6 

3.4E-5 

686 

1.6E-5 

772 

3.0E-5 

858 

6.1E-5 

TABLE  17  MINIMUM  CREEP  RATES  AT  1200°C 


Minimum  Creep  Rate  (s*^) 


Stress  (MPa) 

INb 

2.5  Nb 

5Nb 

9.5  Nb 

3.3  Re 

9.5Nb-3Re 

85 

8.0E-7 

4.2E-7 

86 

6.1E-7 

2.8E-7 

160 

1.9E-6 

169 

1.9E-6 

1.4E-6 

171 

1.4E-6 

8.6E-7 

l.lE-6 

240 

3.8E-6 

254 

4.4E-6 

4.1E-6 

2.6E-6 

257 

4.8E-6 

1.3E-6 

321 

7.9E-6 

338 

8.9E-6 

1.9E-5 

343 

1.3E-5 

1.5E-6 

4.5E-6 

401 

1.7E-5 

423 

1.8E-5 

429 

3.8E-5 

3.0E-6 

8.0E-6 

481 

4.3E-5 

507 

3.5E-5 

514 

1.2E-4 

4.9E-6 

1.5E-5 

592 

7.0E-5 

600 

8.3E-6 

3.4E-5 

686 

1.6E-5 

772 

3.0E-5 

858 

6.1E-5 
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TABLE  19  COMPRESSIVE  CREEP  DUCTILITY  RESULTS 


Creei 

Ductilil 

-y(%) 

Composition 

1 

Creep  Region  of 
interest 

1 

1000°C 

1 

1200°C 

1 

1400°C 

■ 

■ 

i+n 

3.9 

8.7 

21.2 

i+n+m 

9.1 

19.6 

30.5 

(MOo.93  Nbo.07)Si2 

i+n 

1.0 

6.5 

6.2 

i+n+m 

1.4 

11.0 

9.1 

(MOo.85  Nbo.l5)Si2 

I+n 

1.6 

8.2 

18.3 

I+n+m 

1.8 

14.1 

23.0 

(MOo.7lNbo.29)Si2 

I+n 

1.5 

7.1 

12.4 

I+n+m 

2.7 

9.5 

16.3 

(Moo.9Reo.i)Si2 

I+n 

1.7 

1.3 

2.2 

I+n+m 

2.9 

6.5 

5.9 

(Moo.6iNbo.29Reo.i)Si2 

I+n 

1.7 

5.8 

48.0 

I+n+m 

2.0 

9.3 

56.7 
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TABLE  20  ISOTHERMAL  OXIDATION  RESULTS  (1400°C) 


1 

iVeight  Change  (mg/cm^) 

Composition 

20  hours 

120  hours 

+0.5* 

+1.2** 

(MOo.97Nbo.03)Si2 

+0.2 

+0.1 

(MOo.93  Nbo.07)Si2 

+0.4 

+0.3 

(MOo.85  Nbo.l5)Si2 

+0.2 

+0.4 

(MOo.7lNbo.29)Si2 

-1.2 

-1.1 

(Moo.9Rco.i)Si2 

+0.3 

+0.4 

(Moo.6iNbo.29Reo.i)Si2 

+0.6 

-0.2 

ISSt^lUOjO 

ivd  onoQVBVd  oon 


Fig.  1  Parabolic  Rate  Constants  for  SIO2  and  AI2O3 
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STRENGTH  (MPa) 


0  100  700  900  1,100  1,300  1,500 

TEST  TEMPERATURE  (°C) 


Fig.  2  Strength  of  MoSlj  (Ref.  4) 


DENSITY  NORMAUZED  STRESS 
FOR  1%  CREEP  IN  300  HRS 
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Fig.  3  Creep  Strength  of  Various  Intermetalllcs 

and  Nickel-Based  Superalloy  PWA1480  (Ref.  6) 
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Minimum  Creep  Rate  of  MoSlj  and  (Mo, Re)  Si2(Ref.  6) 


100X 


"A"  Phase:  33.2  Mo  -  58.1  Si  -  8.6  Al 
"B"  Phase:  33.0  Mo  -  63.8  Si  -  3.2  Al 
"C"  Phase:  60.9  Mo  -  38.6  Si  -  0.4  Al 
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Figure  6  Backscatter  Electron  Image  and  Electron  Probe  Microanalysis  of  MoSi2  alloy 
with  a  nominal  composition  of  33.3  Mo  -  61.7  Si  -  5.0  Al  (all  in  atomic  %) 
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"X"  Phase:  32.0  Mo  -  34.7  Si  -  33.3  Al 


"Y"  Phase:  18.8  Mo  -  4.8  Si  -  76.4  Al 
"Z"  Phase:  26.2  Mo  -  2.6  Si  -  71.2  Al 


Backscatter  Electron  Image  and  Electron  Probe  Microanalysis  of  MoSi2  alloy 
with  a  nominal  composition  of  33.3  Mo  -  33.4  Si  -  33.3  Al  (all  in  atomic  %) 


"A"  Phase:  32.5  Mo  -  25.9  Si  -  41.6  Al 


"B"  Phase:  60.1  Mo  -  31 .6  Si  -  8.3  Al 
”C"  Phase:  26.2  Mo  - 1.7  Si  -  72.1  Al 


Backscatter  Electron  Image  and  Electron  Probe  Microanalysis  of  MoSi2  alloy 
with  a  nominal  composition  of  33.3  Mo  -  26.2  Si  -  40.5  Al  (all  in  atomic  %) 


"X"  Phase:  32.9  Mo  -  67.1  Si  -0.0  B 
■Y"  Phase:  32.9  Mo  -  0.0  Si  -  67.1  B 


"Z"  Phase:  50.7  Mo  -  0.0  Si  -  49.3  B 


"X"  Phase:  32.6  Mo  -  36,7  Si  -  30,7  Ge 
"Y"  Phase:  34.7  Mo  -  34.9  Si  -  30.4  Ge 
"Z"  Phase:  Epoxy  Mounting  Media 


Backscatter  Electron  Image  and  Electron  Probe  Microanalysis  of  MoSi2  alloy 
with  a  nominal  composition  of  33.3  Mo  -  33.4  Si  -  33.3  Ge  (all  in  atomic  %) 


"A"  Phase:  p.O  Mo  -  39.2  Hf  -  60.8  Si 
"B"  Phase:  33.1  Mo  -  0.05  Hf  -  66.9  Si 
"C"  Phase:  0.2  Mo  -  57.4  Hf  -  42.4  Si 


Figure  1 1 


Backscatter  Electron  Image  and  Electron  Probe  Microanalysis  of  MoSi2  alloy 
with  a  nominal  composition  of  16.7  Mo  - 16.6  Hf  -  66.7  Si  (all  in  atomic  %) 
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"a"  Phase:  26.5  Mo  -  6.7  Nb  -  66.8  Si 


"b"  Phase:  31.5  Mo  - 1.5  Nb  -  67.0  Si 
”c"  Phase:  26.0  Mo  -  6.8  Nb  -  67.2  Si 


"d"  Phase:  4.0  Mo  -  94.5  Nb  - 1.5  Si 


Figure  12  Backscatter  Electron  Image  and  Electron  Probe  Microanalysis  of  MoSi2  alloy 
with  a  nominal  composition  of  27.0  Mo  -6.3  Nb  -  66.7  Si(all  in  atomic  %) 
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100X 


"X"  Phase:  13.1  Mo  - 18.7  Re  -  68.2  Si 


"Y"  Phase:  20.8  Mo  - 11.1  Re  -  68.0  Si 


Figure  14  Backscatter  Electron  Image  and  Electron  Probe  Microanalysis  of  MoSi2  alloy 
with  a  nominal  composition  of  16.7  Mo  - 16.6  Re  -  66.7  Si  (all  in  atomic  %) 
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Fig.  6  Microstructure  of  MoSi2  Prepared  by  Powder  Metallurgy 


Fig.  16  COMPRESSIVE  YIELD  STRESS 
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Fig.  18  1200^0  COMPRESSIVE  CREEP 
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Fig.  21  1000®C  COMPRESSIVE  CREEP 
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Fig.  23  1400®C  COMPRESSIVE  CREEP 
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